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, ^ Abstract 
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Pair creation in a uniform classical electromagnetic field (Schwinger mechanism) 
^ is studied focusing on the time evolution of the distribution of created particles. The 

time evolution of the distribution in time-dependent fields is also presented as well as 
'"^ effects of back reaction. Motivated by the Glasma flux tube, which may be formed 

at the initial stage of heavy-ion collisions, we investigate effects of a magnetic field 
parallel to an electric field, and find that the magnetic field makes the evolution of 
Q a fermion system faster. 
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1 Introduction 

Pair creation in an external field has a long history pioneered by the discovery of the Klein 
paradox [1], which was further investigated by Sauter [2j. They found that tunneling from 
a positive frequency state to a negative frequency state can occur in relativistic quantum 
mechanics. This is a paradox in view of single-particle theory. If one treats it as a problem 
in quantum field theory, the paradox is resolved and particle pair creation from vacuum 
is concluded. 

Another approach to pair creation in an external field is based on Heisenberg-Euler 
effective Lagrangian [3l H] . Heisenberg and Euler found that the 1-loop effective action 
of QED in a classical electric field has an imaginary part, which means a vacuum in an 
electric field is unstable against creation of particle-antiparticle pairs. In 1951, Schwinger 
fully formalized it using the proper time method [S], and thus particle creation in a 
classical electric field is called Schwinger mechanism. He derived the famous formula for 
pair creation probability per unit volume and time: 



where E is an electric field strength and m is mass of electron. 

This non-perturbative particle creation mechanism in strong external fields has a wide 
range of applications; not only original QED problems but also pair creation in non- 
Abehan electromagnetic fields [El El El [IQj and gravitational backgrounds [HI [121 [131 E] • 
Furthermore, the Mott transition (metal-insulator transition) can be discussed in analogue 
with the Schwinger mechanism [151 [IB] ■ 

For e~^e~ pair creation, an electric field strength of the order or above the critical value 
Ec = m?/e ~ 10^^ V/cm is needed, which is beyond current technological capabilities. 
The prospects to observe the Schwinger mechanism at future X-ray laser facilities are 
investigated in Ref.^/Tj. Furthermore, the Coulomb field of superheavy nuclei may cause 
non-perturbative e^e~ creation [181 [IH] • Also the possibility that a supercritical electric 
field is generated at the surface of bare strange stars is proposed [20]. 

In contrast, sufficiently strong field strengths are expected to be obtained in color 
electric fields which are formed just after collisions between high energy particles, such as 
e+e~ and heavy ions. For example, in the framework of Color Glass Condensate (CGC), 
an initial field strength gE ~ ~ 1 GeV^ {Qs is a saturation scale of colliding nuclei and 
(7 is a strong coupling) is realized in heavy-ion collisions at RHIC (e.g. [21]). Thus, the 
Schwinger mechanism is supposed to be a relevant mechanism for hadronic production. 





from a vacuum persistence probability 
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Pair creation by the Schwinger mechanism in e~^e~ colhsions was first studied by Casher et 
aL [22] with the color flux tube model, in which a color electric fleld between qq is treated 
as a constant and classical electric fleld. A considerable number of studies on pair creation 
in flux tubes have followed shifting its stage from e~^e~ collisions to heavy-ion collisions; 
the Lund model [23], flnite-time corrections [21], flnite-size corrections for longitudinal 
direction (the direction of the electric fleld) [23 [20] and for transverse direction [271 EHl EH] , 
perturbation around a background fleld [3U] and effects of back reaction [211 1321 ESI El]- 
Furthermore, kinetic equations incorporating pair creation from backgrounds have been 
investigated based on semi-classical theory [351 ES EI] and quantum fleld theory [381 ES 
110111111121113]. 

The pair creation probability has been applied to a number of studies on particle 
production in high energy collisions. From Eq. ( |1.2[ ) it is recognized that the pair creation 
probability expresses that for collapse of vacuum. Therefore, the pair creation probability 
has no direct information about number of created particle^ and becomes invalid after 
vacuum collapses, or particles are created. Furthermore, there is a possibility that pair 
annihilation occurs. In view of these circumstances, it is more straightforward to treat 
a particle mean number (0, in|aj,ap|0, in), which has direct information about particle 
number, than to treat the pair creation probability. Thus, we deal with mostly a particle 
mean number, or a distribution function in this paper. 

To get a particle mean number, we must quantize charged particle flelds and deflne a 
particle picture speciflcally. In other words, we need to deflne creation and annihilation 
operators of a particle and an antiparticle. For this purpose, canonical quantization is 
more suitable than path-integral quantization. Formulating the Schwinger mechanism in 
terms of canonical quantization has been done by Nikishov [HI HHl |16] as the resolution 
of the Klein paradox. However, how the distribution evolves in time has been unknown 
in that method because only asymptotic states, in and out, are deflned. Also longitudi- 
nal kinetic momentum distributions of created particles have not been obtained because 
kinetic momentum is inseparably connected with time in an electric fleld. 

This obscurity of a longitudinal momentum distribution is closely related to the fol- 
lowing paradox [35]. Suppose that observer A watches pair creation in a homogeneous 
and steady electric fleld, and gets a phase space distribution function of created parti- 
cles n(p), which is independent of the space-time coordinates because the electric fleld 
is constant in space and time. Furthermore, let observer B move relatively to A along 
the direction of the electric fleld. Then, observer B obtains the same distribution with A, 
namely n{p), because a phase space distribution is scalar with respect to Lorentz boost. 
If rewritten in terms of momentum in 5's coordinates p', this distribution takes a differ- 
ent functional form from A^s ; n{p) = n'{p'), except the case that n{p) is independent 
of the momentum in the direction of the electric fleld. For example, if A observes that 
a particle is created with momentum, then B observes this particle is created with 
non-zero momentum. This inequivalence of two observers contradicts the fact that the 
electric fleld is boost-invariant along its direction. That is to say, the uniform electric 
fleld takes the same strength and the same conflguration in both A^s and i?'s frame, and 
thus A and B are on an equal footing. Therefore, the distribution which depends on the 

and dN /dtdx^ are confused in some papers, however, they are different quantity, w can be equated 
with dN/dtdx^ only if w < 1 (see [i51li5] V 
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longitudinal momentum conflicts with the boost invariance of the system. This paradox 
has made a longitudinal momentum distribution in a homogeneous and constant electric 
field unidentified. 

Another obstacle to revealing a momentum distribution in an electric field is ambiguity 
of particle picture. In an external field which breaks Poincare invariance, we lose the 
criterion to define "what is particle" , while in free Minkowski space we have the criterion 
of positive and negative frequency with regard to Poincare invariance. Because of this 
difficulty, longitudinal momentum distributions have been obtained only after an electric 
field is switched off gZl SHI |19] . 

In this paper, we reveal a longitudinal momentum distribution and its time evolution 
by introducing instantaneous positive and negative frequency solutions. We decompose 
the field operator into positive and negative frequency parts instantaneously. This proce- 
dure extracts the information of a particle picture from the field operator. Creation and 
annihilation operators of a particle and an antiparticle are obtained by this decomposi- 
tion. Because the field operator evolves non-trivially under the influence of an electric 
field, the creation and annihilation operators have time-dependence; ap{t),bp{t). This 
time-dependent particle picture enables us to see the time evolution of the momentum 
distribution (0, in|aj3(t)ap(t)|0, in). 

To get a dynamical view of how an intense electric field decays into particles is im- 
portant to gain an understanding of the initial stage of heavy-ion collisions. How quarks 
and gluons thermalize after heavy-ion collisions, especially short thermalization time < 1 
fm/c expected from hydrodynamic calculations (e.g.[HD]), is not fully understood, al- 
though there has been numerous attempts to understand the thermalization process (e.g. 
[5T] and references therein). The CGC predicts that there are strong longitudinal color 
electric and magnetic fields just after collisions of heavy nuclei, the state of which is called 
Glasma [52]. Thus, the color flux tube model is supported by the CGC and attracts re- 
newed interest. In this circumstances, decay of a color electric field due to the Schwinger 
mechanism |35i |36j plays a pivotal role at the initial stage of heavy-ion collisions. Re- 
cently, a thermalization scenario based on the analogy between the Schwinger mechanism 
and the Hawing- Unruh effect has been proposed [53l |5l]. However, as already noted, 
there remains a dearth of understanding of the Schwinger mechanism in spite of its long 



history. In particular, we need to study effects of a longitudinal magnetic fielc *^ , which is 
a new point in the Glasma flux tube. Thus, in this paper we re-investigate the Schwinger 
mechanism as well as effects of a longitudinal magnetic field strictly based on quantum 
field theory, not on the semi-classical approximation |22l l53l [55] or the kinetic approach 
|35l[36lEZ]. 



The remainder of this paper is organized as follows. In the next section, we review 
the formulation of the Schwinger mechanism in terms of canonical quantization for both 
bosons (scalar QED) and fermions (spinor QED) based on mainly Ref . |44] . A kinetic 
momentum distribution in the direction of the electric field cannot be obtained by this 
method because only asymptotic states at t — >■ ±oo, which are labeled by not kinetic 
momentum but canonical momentum, are defined. 



*^The Schwinger mechanism in a magnetic field has been already studied in some papers [131 HH |371 
im [S5] . However, its efi'ects on a system with back reaction have not been studied. 
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In SecjSj we fully reveal the momentum distribution of created particles in a homoge- 
neous and steady electric field by defining a particle picture instantaneously. The resultant 
distributions are static and do not give us dynamical view of pair creation, because the 
background electric field is static. We explain how the longitudinal momentum depen- 
dence of the distribution is compatible with the electric field which is boost-invariant to 
the longitudinal direction. We also show that a magnetic field parallel to the electric field 
suppresses pair creation of bosons and enhances that of fermions. 

In Secjlj we deal with non-steady electric fields, strength of which is at first time and 
non-zero later. This set-up of electric fields enables us to investigate the time evolution 
of distributions. The particle picture associated with instantaneous positive and negative 
frequency solutions offers us physically reasonable behavior of the distributions, although 
it has been criticized in cosmological particle creation. Also the time evolution of total 
particle number density and charge current is shown. 

In SecjS} we consider back reaction, which is an electric field generated by particles 
created from an initial electric field. The distributions show plasma oscillation, which is a 
classical dynamics. They are also influenced by quantum effect, namely Bose enhancement 
or Pauli blocking. We estimate the time when the initial electric field decays into particles. 

Sec|6] is devoted to summary and discussion. 

In appendix, 2-component formalism of the Klein-Gordon field and a method to solve 
the Dirac equation in an electromagnetic field are reviewed. 

Throughout this paper, we use the units in which h = 1, c = 1 and fc^ = 1. 

2 Canonical quantization in an external electromag- 
netic field 

In this section, we review canonical quantization in a constant electromagnetic field, and 
derive a particle mean number [H] (see also O E]), however a longitudinal momentum 
distribution of created particles cannot be obtained in this section (longitudinal means 
the direction of the electric field) . 

2.1 Bosons in a constant electromagnetic field 
2.1.1 Pure electric field 

We treat charged scalar fields interacting with a classical background electric field: 

C = {df,- ieA^)<l)^ ■ {d^" + ieAf'^ - m^(j)^(f). (2.1) 

At first, we take a constant and homogeneous electric field directing along 2;-axis, and 
choose the gauge potential as 

Af" = {0,0,0, -Et). (2.2) 

This gauge greatly simplifies our calculation because it is independent of the space co- 
ordinates*^ and thus we use the gauge A^ = throughout this paper. The equation of 



*3 



Other choice of gauge, e.g. — {—Ez, 0,0,0), brings some difficulties (see [571 [55] 1. 
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motion for is 

[d^ -d^- (^3 + ieEty + m'](P = 0, 
where d'^ = d1 + dy. The field may be expanded as 



(/.(t,x) 



(2vr; 



3/2 



(2.3) 



(2.4) 



because the equation of motion (2.3) does not depend on the space coordinates exphcitly. 
The mode function 0p(t) obeys the equation 



+ {eEt + p^Y + mf^ 



in which transverse mass mx is introduced as 



= + = + pi + Py. 

Note that the role of transverse degrees of freedom is only shifting the mass. 
Changing variable as 

^ = J^{eEt + p,) 



converts Eq.(2.3) into a Schrodinger-like equation 

' d^ 1 



0p(t) = 0, 



in which a is dimension-less parameter defined by 



2^' 



If we choose two solutions of Eq.(2.8) ^^{t) satisfying conditions 

z^;(t)^+0p(t) = i 



z_0;(t)-_0p(t) 



dt 
d 
dt 



z+0;(t)-_0p(t) = o, 



(2.5) 



(2.6) 



(2.7) 



(2i 



(2.9) 



(2.10) 



then the mode functions defined by ±fp{x) = -ti^p(t)e*P'''/(27r)'^/^ fulfill the orthonormal 
conditions 



;+/p,+/q)=^'(p-q) 
-/p,-/q) = -5'(p-q) 
(+/p,_/q) = 0, 



(2.11) 
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where the inner product (/, g) is defined by 



if,g) = t d'xfdog. 



(2.12) 



Because of the orthonormal conditions (2.11), the field operator may be decomposed as 

(2.13) 



-ap + -(f)p(t)- 



ip-x 



(2vr)3/2 P -^P^ '(27r)3/2"-PJ 
j d^p[+fpix)ap + Jp{x)blp], 



hi 



(2.14) 



and relations 

-p -(-/p,0) 

hold in the same way as the free field. Imposing the canonical commutation relation 
[0(t, x), 7r(t, x')] = ^^^(x — x'), where 7r(x) is canonical conjugate momentum 7r(x) = (p^x) 



and using Eq.(2.14), one yields commutation relations 

[op, a]j] = [hp, 6y = ^^(p - q), others = 0. 

The charge operator Q and the canonical momentum operator are expressed as 



(2.15) 



Q = ie d X 



d 

di' 



(2.16) 



e j d^p^a^ap - h^hp) 



d^p p\alap + hl,hp) (z = 1,2,3). 



(2.17) 



In the derivation of these equations, we have used Eq.(2.11) and subtracted the infinite 



vacuum contribution by normal ordering. Because of Eqs.(2.15),(2.16) and (2.17), aj, and 



Op can be interpreted as the creation and annihilation operator of a particle with charge 
+e and canonical momentum p respectively, and and hp as the creation and annihila- 
tion operator of an antiparticle with charge — e and canonical momentum p respectively. 
Notice that Q can be diagonalized by the creation and annihilation operators owing to the 



orthonormal condition (2.11), and P* is diagonalized owing to the orthonormal condition 



(2.11) and spatial translational symmetry as well. These are independent of an actual 



choice of the solutions of Eq.(2.8) i(/>p(t). 



One set of solutions of Eq.(2.8) satisfying the orthonormal conditions (2.10) can be 



expressed in terms of a parabolic cylinder function D^{x) |59l |60] as 



{2eE)t 



r/^_.,_|(ef^O 

5 2 

rD* . i(ef'0- 



(2.18) 
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Choice of these solutions is not unique. There remains the freedom of the Bogohubov 
transformation: 

+Mt) = "p4^pW + P*p-Mt) ^2 19) 

_0p(t) =a;_0p(t)+/?p+0p(t) 

with the constraint |app — |/3pp = 1. Since these Bogohubov transformed solutions -t^^p 



also satisfy the conditions (2.10), quantization procedure can be done choosing any -y/)p 



from these infinite sets of solutions, and thus there are infinite degrees of freedom how to 



expand the field operator (2.13), from which the creation and annihilation operators are 
defined. Accordingly vacua defined by ap|0) = &p|0) = exist in infinite way and any two 
of them are not equivalent [611 E2] . 

The identical arbitrariness also exists in quantization of the free field. But in that 
case, it is natural to choose mode functions as eigenfunctions of id^, that are positive and 
negative frequency solutions with regard to Poincare invariance. 

In an external field, decomposition of mode functions into positive and negative fre- 
quency solutions is not straightforward because the Lagrangian depends on space-time 
coordinates explicitly. In other words, particle concept becomes ambiguous due to inter- 
action with the external field. 

Therefore, we need a criterion to define "particle". Now, we employ "asymptotic 
WKB criterion" [HI [7|, which leads to the same result as obtained in the proper time 
method by Schwinger. In this criterion, we determine a particle picture with the same 
idea as in a scattering problem in field theory. Namely, decide a particle picture from 
the behavior of mode functions in asymptotic regions, infinite past and future. But now, 
there is no asymptotic free region and the mode functions do not approach the plane wave 
solution exp(=Fipx), because we treat the steady electric field which lies from infinite past 
to infinite futur^*^ Thus, we define positive and negative frequency mode functions in 
such a way that the positive frequency mode function approaches exp[iSci(t)] instead of 
exp{—ipx) and the negative frequency one exp[—iSci(t)] instead of exp(ipx) in asymptotic 
region*^, where Sd(t) is a classical action of a charged particle under the electric field 
A^" = (0,0,0,-Et): 



-\eEt^ - a In (^y/eEt^ {t +oo 
\eEt^ + a\n (yeE\t\^ {t -oo 



(2.20) 



One can confirm that the mode functions (2.18) behave like as 



4^p(t) ^ {y/eE\t\)-'2 exp{-|eEt2 - mln(v^|t|)} 

^2.21^ 

_0p(t) ^ {VeE\t\)~h exp{{eEt^ + zaln(v^|t|)} 

**Instead of using the asymptotic WKB criterion, one can employ an electric field which asymptotically 
vanishes at infinite past and future, and obtain the same result with us |45l 146 j . 
*^Notice that a classical action of a free particle is —px ~ — -^/p^ + vnP-t + p • x. 
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at t — > +00 modulo a constant factor, using the formula 



z^'e (I argz| < |7r) 



D,{z) ^ ^ \ z^e--^'^ - ^f'^'z-^-^el^'" (^tt < |argz| < \r) (2.22) 



|^|— »oo 



^'^e"!^' - ^^e-'^^^z-'^-^e^" (-fvr < | arg^] < . 



Therefore +(j)p{t) and -4'p{t) given by Eq.(2.18) are settled as the positive and negative 
frequency solution, respectively. 

At t — > —00, however, these functions do not behave like as :i4>p(t) ~ exp[±iSci(t)] 
and another mode functions do. This is the very consequence of pair creation; particle 
picture in infinite past and that in infinite future are not equivalent. Thus, let ±4>p{t) 



defined by Eq.(2.18) be ±0p"*(t). On the other hand, one can find that 



^ (2eE)3 2 

_C(t) = ^D_,,_i(-ef^O 

P ^ ^ (2e_B)3 2 ^ ''^ 



(2.23) 



behave like as i(/'p(t) ~ exp[±iSd{t)], respectively, at t ^ —00. 

/,out/+\ I a* j.out/ 



±0p and ±0°"* are related by the Bogoliubov translation as 



(2.24) 



with coefficients 



r(i -m) 

One can confirm the equation 

|ap|'-|/5p|' = l (2.26) 

which guarantees the unitarity, holds. 

We have chosen two sets of mode functions ±0p and ±0°"* from infinite candidates 
using the asymptotic WKB criterion. Then, we can define two distinct particle pictures 
by expanding the field operator using ±fp^°^^{x) = -t0p''°"*(t)e*P'''/(27r)^/^: 



(2.27) 



0(x) = J dMj'^ix)a^ + j'^ixW^i] 
and two distinct vacua are introduced as 



a!,"|0,in) = 6i;,'|0,in) = (2.28) 



"p I"! "p 

a°"*|0,out) = 6°^*|0,out) = 0. (2.29) 
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Using (2.11), (2.14), (2.24) and (2.27), the relations between creation or annihilation 
operators of in and out are derived: 



^out 



out 

p 



(,outt 

p 



-jp 



out 



(2.30) 



dMi-fT,J'^)a'^ + UT,-Ob 



Now we can calculate a mean number of particle pair produced by the electric field 



using the relations (2.30). Assume now that the state is |0, in), which is the vacuum of the 



particle picture defined at t — > —oo and which is invariable because we use the Heisenberg 
picture. By definition (2.28), the mean number of in-particles is 0: (0, in|ap"''ap |0, in) = 



(0, in|6p"l'6p |0, in) = 0. However, we must consider that particles observed really are out- 
particles, and they are condensed in the in-vacuum. 

We define the distribution function of particles and that of antiparticles by 



Up = (0,in|ap 



outf out 
Up 



0,in) 



and 



Hp = (0,in|6°"*t5°"*|0,in) 



V 
V 



(2.31) 



(2.32) 



respectively. Because aj,ap (&p&p) is the operator which represents the total (anti)particle 
number in whole space, their expectation values diverge proportionally to space volume. 
Therefore, we have divided them by infinite space volume V = jd^x = (27r)^(5^(0) and 
converted them into the (anti)particle number per unit volume, which is finite. Further- 
more (27r)^, which is really {2nh)^ and represents the unit of phase space volume, has been 
multiplied to get the (anti)particle number per the unit of phase space volume (27r/i)^, 
that is a phase space distribution function. 



Using Eq.(2.25) and (2.30), the distribution functions (2.31) and (2.32) become 



-2-Ka 



1/5- 



= exp 

|2 

-pl i 



TTrrin 



eE 



(2.33) 



where 



holds for general ctp and (3p because of the charge and the momentum 



conservation. Thus, here and henceforth we call rip particle pair distribution function, 
which means the number of particles with momentum p and that of antiparticles with 
momentum — p simultaneously. 

Note that Up is independent of the longitudinal momentum p^. In view of Lorentz 
invariance, this fact is reasonable since the system with an electric field directed along 
2;-axis is invariant under Lorentz boost to z-direction. However, it is strange that there 
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equally exist particles with every momentum including opposite direction to the electric 
field. 

This paradox comes from confusion of canonical momentum with kinetic momentum. 
Now the subscript p of a-p and hp 



represents canonical momentum because Op and hp 
diagonalize the canonical momentum operator (2.17). Because a°"* and 6°*^* have been 

pout;' 
'P 

+CXD, ttp"* and 6°"* are the operator which an- 



introduced by (2.27) in which all time-dependence was taken by ±/p"*(a;) and they were 



determined to approach exp[±iS'ci(t)] at t 
nihilate a particle or an antiparticle with canonical momentum p at t — > +oo, respectively. 
With the gauge = (0, 0, 0, —Et), canonical momentum p, which is a gauge-dependent 
quantity, has relation to kinetic momentum 11 as 

= n, - eEt, (2.34) 

which can be explained that kinetic momentum increases linearly with time while canon- 
ical momentum is a constant of motion with regard to translational symmetry. (For 
transverse direction, canonical and kinetic momentum are identical because there is no 
external field in that direction.) If we take a limit t — > +oo for fixed Pz, kinetic momentum 
Hz approaches infinity. Thus, for any value of pz, rip represents the number of particles at 
t — *• -|-oo, when kinetic momentum of particles becomes infinity due to acceleration by the 
electric field. Of course in reality, there are particles with finite kinetic momentum at any 
time since pair creation happens constantly, however, they are invisible by this method 
because we took the limit t — > +oo for fixed Pz- If one want to see a kinetic momentum 
distribution of created particles, a particle picture must be defined at each time not only 
asymptotic region. This will be done in SecjSj 

Although the z-component of canonical momentum is not a physical value, we can 
interpret it as an indicator of the time when pair creation occurs. Letting the time of pair 
creation be to and assuming kinetic momentum of the created particle is at this time, 
Eq.(2.34) indicates to = —Pz/eE. Because the system is static and pair creation happens 
constantly, there is no special value for to- This is the reason why Up is independent of 

Pz- 



Also the pair creation probability can be derived by our method. Although we have 
little interest in the pair creation probability as a physical value, we show its derivation 
in order to clarify meanings of each calculation. 



From Eq.(2.30) one can find the relation between |0,in) and |0,out): 

^3 



|0,in) = JjA^p-^/^gxp 



(27r)3 (3p 



|0, out) 



(2.35) 



1 



l-^pl = I«PI 

The in-vacuum is a two-mode squeezed state of the out-particle pairs. Using this relation, 
we can calculate a vacuum persistence probability as follows: 

V 



Po = |(0,out|0,in)|2 = Y[ lATpl-i = exp 



(27r)= 



(fpln (l 



— 2na 



(2.36) 



The last expression's index diverges in two ways. Because of this divergence, (0, out|0, in) 
is zero, which means that the Fock space built on |0,in) and that on |0,out) are inequiv- 
alent 1621. 
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The first divergence occurs because of infinite volume: V = (27r)^5^(0). That is to 
say, because the pair creation probability per unit volume is not 0, pair creation happens 
without fail in infinite volume system. 

The second one comes from J dpz, which corresponds to the fact that the time interval 
when pair creation can occur is — oo to +00, because pz is the indicator of the time 
when pair creation happens. The relation dp^ = eE dto = eET follows from 
to = —pz/eE, where T expresses whole time interval, which is infinite. 



We can now rewrite the vacuum persistence probability (|2.36|) as 

VTeE 



= exp 



jd^PT ln(l + e-2™) 



(2.37) 



from which the pair creation probability per unit volume and time can be read as 

eE — ^-Ti^+i 



w 



(27r)3 

{eEf 
f27r)3 



d pt 



^ n 



-2Txna 



n=0 



\n+l 



E( 1 ) Trnm 




(2.38) 



n=0 



2.1.2 Collinear electric and magnetic field 

We take a collinear electric and magnetic field E = (0, 0, £"), B = (0, 0, B) motivated by 
the Glasma flux tube |52]. We use the gauge potential = (0, —By, 0, —Et). 
The Klein-Gordon equation is 

[dl - {dl + ieByf - dl - (^3 + ieEtf + m2]0(t, x) = 0. (2.39) 

This can be solved by separation of variable: (l){t,x) = '^{t)x{y)i^G^^^"'^^^''^'^ ■ Then, ipit) 
and xiy) obey equations 

— -{eBy + p,f + C 



dy^ 



d^ 
Jt 



2 +{eEt+Pzf + + C 



ifiy) = 

x{t) = 0. 



These can be rewritten as 



d^_rf J_ 

h — H ^ 

de 4: 2eE 



x{t) = 



with variables 



2 / 2 

— {eEt + p,), y] = \J —{eBy + p^ 



(2.40) 
(2.41) 

(2.42) 
(2.43) 

(2.44) 



Eq.(2.42) has bounded solutions only if C, takes discrete value, = (2^ + l)e-B (n 

T feB\'^ 1 



0, 1, • • ■ ), and a solution is 



27r V TT / 



(2.45) 



12 



which is labeled by the Landau level n and normalized as 

L 



dy(j)n{y)(l)m{y) 



where L is the linear size of the system: V = L^. 



(2.46) 



Eq.(2.43) takes the same form as that in the pure electric field case (2.8) with the 
replacement m^ + (^„ — > rn^+p^. Thus, we can get the solution of (2.43) from the solutions 
in the pure electric case, (2.23) and (2.18), by the replacement a — ™ ^^'^ ^ " — 



2eE 



2eE 



{2eE)^ 



[-en) 
[-em 



or 




{2eE)^ 



[em 



(2eE)^ 



(2.47) 



where p is the abbreviated expression for {px,pz,n). 

Consequently, we can connect the pure electric field case and the coUinear electromag- 
netic field case by the replacements 



V2i^ 

jdpy 



27r 



m^+(2?i+l)eB 



(2.48) 



2eE 2eE 

The pair distribution function in the presence of the collinear electromagnetic field is 

7r[m2 + (2n + l)e5] 



rip = e 



-2'Ka'n 



exp 



eE 



(2.49) 



from which one sees that a longitudinal magnetic field makes particles "heavy" and sup- 
presses pair creation. 



We can also derive the pair creation probability by the replacements (2.48). The 



vacuum persistence probability (|2.36|) changes to 

L2 



Po = |(0,out|0,in)p = exp 



(2vr)^ 



jdp, jdpjn{l + e-^-^-) 



(2.50) 



The divergence from J dp^ is treated in the same way as that in the pure electric field 
case: / dpz = eET. 

New divergence comes from J dp^ because the integrand is independent of px in this 
case. This is related with degeneracy of the Landau orbit. Comparing Eq.(2.40) and the 
Schrodinger equation of one-dimensional harmonic oscillator 



J_d^ 
2m dx"^ 



-muj'^x^ -|- e ) = 0, 



(2.51) 



we can notice that ?/o = is something like a center of circular motion. Because 

a center of circular motion does not influence physical result, the distribution function 
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(2.49) and the integrand of (2.50) is independent of Px- Using the relation = — the 



divergence from f dp^ can be treated as / dp^ = f d^—eByo) = eBL. 
Then, the vacuum persistence probabihty is rewritten as 



= exp 



^5^1n(l + e-2™") 



(27r)^ 



(2.52) 



from which the pair creation probabihty per unit volume and time can be read as 



w 



(27r)2 ^ ^ / ' 
e^EB ^ (-1 



71=0 1 = 
00 , 



(2.53) 



(2vr) 



1=0 



I 



-e 



2sinh7r/|' 



which is also suppressed by the magnetic field. This result was first derived by Popov 
with the imaginary time method |55i|. 

2.2 Fermions in a constant electromagnetic field 

Next, we study pair creation of fermions which obey the Lagrangian 

C = ^P [ij^{d^ + leA^) - m] i), (2.54) 
in the same way as that of bosons. 

2.2.1 Pure electric field 

The equation of motion derived from the Lagrangian is 

Yl^{%d^ - eA^) - m] ^{x) = 0, (2.55) 
of which solutions under the gauge potential = (0, 0, 0, —Et) are ±-?/'p^(a;) and -iip'^^{x) 



given by Eq.(B.28) and (B.29 



is the same as that in Sec, 2.1 



respectively (see Appendix [B]). The meaning of in and out 
The subscripts t, | represent spin directions. 



These solutions are related by the Bogoliubov transformation: 



+V;i^,(x) = ap#°f (x) - (3;-%fix) 
_<(x) = al^r^fix) + /3p+^°f (a;) 

where the coefficients are 

\/27ra^ 



=T,i), 



(2.56) 



r(l -ia) 

These coefficients satisfy the relation 



(2.57) 



(2.58) 
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Notice that the sign on the left-hand side is opposite to that of boson [Eq.( |2.26 )] as a 
consequence of spin statistics. Because a distribution function is given by |/5pp not only 



in the boson case but also fermion [this will be shown later], Eq.(2.58) guarantees Pauli's 
exclusion principal: rip = |/5pp < 1. 



The solutions iV'™ and i^"" satisfy orthonormal conditions 



where the inner product is defined as 



as 



as 
ps' 



h^qs' 







as = in, out 



(2.59) 



(2.60) 



If one expands the field operator in two ways: 



^(x) = J2 d'p Upsi^Ks + -^ii^w 



inf 



(2.61) 



and imposes canonical anti-commutation relation {^/'(t, x), 7r(t, x')} = ^^^(x — x'), (7r(x) 
itfj^x)), following equations 



Ks. <2} = {&ps, &qi} = '^'(P - q)'^..' , others = 



=T,i 



P'=Y. I P'Ks<s + &p>pJ = 1, 2, 3) 



(2.62) 
(2.63) 

(2.64) 



(as = in, out) are obtained with the help of the orthonormal conditions (2.59). With these 
equations, we can interpret a™^ [6qj'] as the operator which creates a particle [antiparticle] 
with canonical momentum p and spin s at infinite past, and flq"*^ [^q"*^] as that at infinite 
future. 

The relations 

[Wps 



ar.* = (+Cf,^) 



Lint 



(2.65) 



OipO_pg Ppfl 



ps 



are followed from Eq.(2.56), (2.60) and (2.61) 



Now, we can calculate a distribution function of created pairs observed at time of +00: 

(2.66) 



np. = (0,in|a°f a°";|0,in)%^ = (0, in|O-p|0' in)-^^""^' 



V 



V 
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Of course, this is independent of spin states because there is no magnetic field. 



Notice that the distribution (2.66) is identical to that of boson (3.17). It seems strange 



that distribution functions do not reflect spin statistics. This paradox can be understood 
as follows: we now see particles with canonical momentum p at infinite future, and as 
mentioned earlier, kinetic momentum of these particles are infinite due to acceleration 
by the electric field at that time. Therefore, spin statistics vanishes because of infinite 
energy. 



Next, we calculate the pair creation probability. Using Eq.(2.65) one can derive a 
relation between |0,in) and |0,out): 



io,in) = nA/; 



-1/2 



pT 



(27r)3 /^p outfTOutt 

^ V a* ''-pT 



UK 



-1/2 



qi 



q „outtjoutt 



iA/pTi = iA/; 



Pi I 



1 



-27ra ' 



This relation yields the vacuum persistence probability: 



Po = |(0,out|0,in)p 



nn 

s=hl P 



Or 



exp 



from which the pair creation probability is obtained as 

2 



w 



(27r)3 

2eE 
(27r)3 



^ n=0 

oo _ 

E- 



2iTa 



-Knm 

■e 



n=0 



Of course, this is identical to the Schwinger's result. 



|0, out) 
(2.67) 



(2.68) 



(2.69) 



2.2.2 Collinear electric and magnetic field 

Similarly to the boson case, we can connect the pure electric field case and the collinear 
electromagnetic field case by the replacements 



Jdpy 



2eE 



27r 
L 



(2.70) 



+2neB 



n — 2eE 
,2 



m^+(2n+2)eB 
2eE 



(n 



0,1, 



where n and m are Landau level. Compared with the boson case (2.48), energy levels are 
split due to the spin-magnetic field interaction. 

Thus, the distribution in the collinear electromagnetic field is 



ri- 



ps 



exp(— 27ra^ 



4,i), 



(2.71) 
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in which p is the abbreviated expression for {px,pz,n). 

Notice that, also in this case, the magnetic field behaves like mass and suppresses pair 
creation. On the other hand, spin-magnetic field interaction makes a particle of |-mode 
lighter. Especially, n = and j-mode is not suppressed by the magnetic field because % 
is independent of B. 



The replacements (2.70) changes the calculation of the vacuum persistence probability 
as follows: 



b = exp 



exp 



(2vr) 



VTe^EB 



(27r) 



-27ra: 



■•) 



,1 ri 



(2.72) 



in which the divergences from J dp^ and J dp^ are treated in the same way as that in 



Sec. 2. 1.2 Then, the pair creation probability is 

e'^EB 



w 



{2ny 
e^EB 



J2 5^1n(l-e-2-<) 



s=T,i n 



(27r)2 2^ 2^2^ f 

(P'EB X ^ 1 TTim^ B 

— — > -e cotliTr/ — . 
27r 2 ^ E 

^ ' 1=0 



(2.73) 



This w increases with increasing B*^ while that of boson [Eq.(2.53)] decreases. This 
is because (i) pair creation of fermion in n = and j-mode is not suppressed by the 
magnetic field, while that of bosons in any Landau level is suppressed exponentially, and 
(ii) the number of modes degenerating in one Landau level per unit area is eB/2'K, which 
is proportional to the magnetic field strength, as is known in quantum dynamics (see for 
example |S3])- 

This phenomenon is common to other bulk quantities: expectation values of the Hamil- 
tonian, the number operator or the current operator, etc. For example, the in-vacuum 
expectation value of the total number operator for fermion out-particles is 

(0,in|iV|0,in) = y"^P.y"^P.(0,m|a°f a-;|0,in) 



--hi 



2n 
eB 



dpz 

27r''P^ 



(2.74) 



V 



dpz 
2t\ 



^This behavior was first pointed out by Nikishov 
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which is approximately proportional to B^J^ From Eq.(2.74) we can confirm that the 
degeneracy of one Landau level per unit area is surely eB/2T[. 



3 Longitudinal momentum distribution 

In the previous section, we have calculated the distribution function of created particle 
pairs using the asymptotic WKB criterion. However, a kinetic momentum distribution in 
the direction of the electric field (longitudinal direction) could not be obtained because 
we saw only the asymptotic state, in which longitudinal momentum becomes infinite due 
to acceleration by the electric field. 

To study a longitudinal momentum distribution in the presence of the constant electric 
field, we must define a particle picture at each instant not only at infinite past or future. 
That is, to expand the field operator by instantaneous positive and negative frequency 
solutions at each time. 

We define instantaneous positive and negative frequency solution at time t = 1^ as the 
plane wave solutions under the pure gauge A = A(to), where A(t) is a gauge potential 
for an electric field E(t) (we use the gauge A° = 0). For a scalar field in the constant 
electric field J^if) = —Et, the instantaneous positive and negative frequency solutions at 
t = to are, respectively. 



(*0)| 



1 



XI 



V(27r)32o;p 



3=Fja;pt+j(p-Eto)-x 



(3.1) 



in which lo.,, 



+ m\ 



*8 



<j){x) 



. Then, the field expansion is 



(3.3) 



where ap(to) and &-p(to) are annihilation operators defined at t = t^. Because the mode 
functions (3.1) are not solutions of the field equation in the electric field, the creation and 



annihilation operators are forced to have time-dependence, which enables us to see the 
time evolution of the momentum distribution. 



*7r 



The integral J'^n-ps in Eq.( |2.74 l diverges because now Ups is independent oipz- Finite expectation 



values will be obtained in Secj4] where an electric field is imposed during finite time. 
*^In Ref.SOl Emilia] 



±/p(a;) 



1 



V(2^)32a;p(t) 



exp 



=Fi / (ii'a;p(t') + ip • X 



(3.2) 



where i^p{t) — -^( p — eA(t))^ + m? is used as instantaneous positive and negative frequency solutions 
instead of Eq.(3.1). Equations (3.11 and (3.2 1 give the essentially same particle picture although the 



meaning of p in (3.1) and ( |3.2| are different; p is kinetic in (3.1) and canonical in ( |3.2[) [One can confirm 
it by calculating the canonical momentum operator]. Differences between us and Ref . |4oi I¥T1 are as 

follows. In Ref . [^] |4T1 |42j |43] , a kinetic equation is investigated based on quantum field theory. However, 
it is practically easier to treat the field equation itself, which is a second-order differential equation, than 
to treat the kinetic equation with a non-Markovian source term, which is derived from the field equation 
and is an integro-differential equation. Furthermore, the employed regularization methods are distinct. 
In our method, the Hamiltonian is regularized only by normal ordering, and the expectation value of the 
current operator is finite without regularization. 
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The idea of instantaneous positive and negative frequency was used to study particle 
creation in expanding universes. However, it encounters some difficulties and is strongly 
criticized [Ml [ll]- The main subject for the criticism is the behavior of a Bogoliubov 
coefficient f3p as a function of p. For the Robertson- Walker model, jSp falls off with 
negative power of |p| for high momentum p. Especially for anisotropic case, jSp behaves 
like jSp ~ |p|~^ for high momentum, and therefore the particle number per unit volume 
/ (f^l'^^pl^ diverges [H]. This is physically unacceptable. 

Because of above-mentioned problem, other approach is believed to be a correct one 
to treat cosmological particle creation. It is adiabatic approach [HSl [131 [II]- In the 
adiabatic approach, a physical Fock space is fixed by expanding the field operator by 
the exact solutions whose asymptotic form for large mass or momentum is a WKB-type 



solution^ Then, one ceases defining particle picture at intermediate time and treats 
mostly field observables like as the energy-momentum tensor and the charge current, 
which are renormalized by adiabatic regularization. For an electric background field, this 
approach has been used by Kluger et al. [321 155] . 

In spite of these circumstances, we apply the idea of instantaneous positive and nega- 
tive frequency and investigate a particle number at each instant of time. That is because, 
for an electromagnetic background, (3p falls off rapidly for high momentum if one take an 
appropriate electric field configuration (see Secjij). [In this section, /3p does not fall off for 
high longitudinal momentum. It is not a defect of our method but a natural result since 
we take a constant electric field.] Thus, the application of the instantaneous positive and 
negative frequency to particle creation in an electromagnetic field can avoid the criticism 
subjected to cosmological particle creation. 

Of course, as noted earlier, a definition of particle in the presence of a pair-creating 
background is rather ambiguous. The use of the particle picture defined by the instanta- 
neous positive and negative frequency solutions is merely an ansatz. An ansatz must be 
justified by physically reasonable results. Actually the ansatz of the instantaneous parti- 
cle picture will provide us with convincing behavior of distribution functions. Especially 
in Secjij we will get distributions which yield a finite total particle number. 

A merit of this particle picture is that the divergent part of the energy-momentum 
tensor can be removed only by normal ordering; no further regularization is needed unlike 
the adiabatic approach. This good point benefits even those who do not trust the concept 
of particle number at each instant and accept only field observables. One can interpret 
defining a particle picture at each time as a means to identify zero-point energy. Fur- 
thermore, energy conservation is automatically and manifestly guaranteed in this particle 
picture because the energy operator at each time takes the same form as that of the free 
field: 

P\t) = J d?p LOp [at (t)ap(t) + bp{t)bl{t)] , (3.4) 

and the Bogoliubov transformation, which relates creation and annihilation operators of 
different times, conserves the commutation relation: bp(t)b^p(t) — : 6p(t)6|,(t) : = const. 
This is not the case generally in the adiabatic approach, in which subtracted part is time- 
dependent. In this and next section, however, total energy of the system is not conserved 
because there is an external electric field which provides the system with energy. If one 



*^This procedure is identical to the asymptotic WKB criterion we have used in Secj2j 
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takes account of back reaction and treat an electric field as a dynamical variable, total 
energy will be conserved (Secjsj). 



3.1 Bosons in a constant electromagnetic field 
3.1.1 Pure electric field 

We re-investigate the distribution function of created particles in the constant electric field 
using the particle picture defined by the instantaneous positive and negative frequency 
solutions. That is to say, we calculate the pair distribution function 



np{t) = (0,in|aT (t)ap(t)|0,in) 



(2vr) 



V 



(0,in|6lp(t)6_p(t)|0,in) 



(2vr) 
V 



(3.5) 



The instantaneous particle picture is introduced and related to the in-particle picture by 
the expansion: 



inf] 



t=to 



(3.6) 



The creation and annihilation operator of in-particle and that of time to are related by 
the Bogoliubov transformation whose coefficients «p,q(t), /3p,q(t) are given by 



Jd^q 
Jd^q 



it), 



«;,qW-/f^(^)+/5p,q(t)+/. 



(3.7) 



However, this expansion is not definite because the coefficients has time-dependence. This 
ambiguity of the expansion comes from the fact that the Klein-Gordon equation is the 
second-order differential equation. Accordingly, to make the expansion determined, we 
put the additional conditions 



-/p(^)l*=to 



Jd^q 
Jd'q 



«P,q(to)+/i*°Hio, x) + (3;^^{to)-fi''\to, x) 

«;,q(i0)-/^'Ht0,x)+/3p,q(t0)+/i*''^(t0,x) 



(3. 



where dot denotes time-derivative. 



Incidentally, these conditions, (3.7) and (3.8), tell us that the particle number defined 



by the instantaneous particle picture at time t is identical to which is obtained if one 
switches off the electric field at time t. This accordance is merely the matter of calculation; 
we do not switch off the field but only define a particle picture at each time. 

In our calculation, it is convenient to use the 2-component form of the Klein-Gordon 



field (see Appendix |A|). If one put the same condition as Eq.(3.7) on the 2-component mode 
function ±Fp\ which is made from -t/p^ by Eq.(A.2), the condition (3.8) is automatically 
satisfied. 

Furthermore, we can get the Bogoliubov coefficients easily using the 2-component 



formalism. Because the inner product expressed by the 2-component form (A. 5) does not 
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contain time-derivative, the Bogohubov coefficients can be extracted from Eq.(3.7) by 



the inner product: ap,q(to) = {+f^°\+fn)^ while the inner product of the 1-component 



Eq.(3.7). 



form (2.12) contains time- derivative and cannot extract the Bogoliubov coefficients from 



The 2-component forms corresponding to ±fp{x) = -t//)™ (t)e*P'^/(27r)^/^ are 



1 / Trnxj)^ + i^^ 



(3.9) 



where ^f/'p (t) are given by Eq.(2.23), and those to the instantaneous positive and negative 
frequency mode functions ±fp°\x) [Eq.(3.1) 



are 

1 I m + uj. 



p 



2m V m-uj, 



p 



m — UJr. 



2m \ m + Up 

Using these expressions, one can calculate the inner products as follows: 



(3.10) 
(3.11) 



2uJr 



e i e 



(^p+p,)DVi(-e^' 




Pz 



+ e^'V2^D' 



(2eE)V4 ^ 

X 5^(p - q - eEto) 
= ap(to)5^(p - q - eEto) 

2m I d'xJ'^''^\to,^)cxs+F'^{to,^) 



^P+C-eEto(^o) - n0?-eEio('^o) ^'(p - Q " eEto) 




Pz 



(2eE)V4 ^ 

X 5^(p - q - eEto) 
-(3;{to)6'\p-q-eEto) 




(+/^"^-/q) =2^ y"rf^X4i^(*°)t(i^^x)a3^^(to,x) =/3p(to)5'(p-q-eEto 
~fp'\-0 = ^rnj rf3a;_F(*°)t(t^,x)a3^'^(to,x) = -a;(to)5'(p - q - eEto), 



(3.12) 
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from which the relations 



+fpi^) — "p+eEto(^o)+/p+eEto(^) + /^p+eEto (^o) -/p+eEio 



{to) 



-fp{x) — ttp+eEto(^o)-/p+eEto(^) + /^P+eEto (^o)+/p+iEto 

are obtained. One can check that the equation 

i«p(^o)r-i/9p(to)i' = i 



{to) 



(3.13) 



(3.14) 



holds using the explicit form of ap(to) and /3p(to), or the orthonormal conditions for each 
mode functions. 



Inserting (3.13) into (3.6), we get 



ap(to) - «p(^o)ap-eEto + /^p(^o)&!^'p+eEto 
blpito) = «;(to)&'L";+eE*o +/^;(^oXeEto- 



(3.15) 



If one calculate the ^-component of the canonical momentum operator using the second 



expression in Eq.(3.6) 



p3 = / (fp [{p, - eEto)al{to)aM + {p, + eEto)bl{to)bp{U 



(3.16) 



is obtained, from which a'^p{to) [bp{to)] is identified as the operator which creates a particle 
[antiparticle] with kinetic momentum p at time to- 



Now, we can calculate the pair distribution function at time t. Substituting (3.15) 



into (3.5) yields 



^pW = l/3pWP 



e 2' 



[UJr 



p.)D_^.{-e^" 




Vz 



-64' 




Pz 



(3.17) 



in which p^ denotes, as noted above, kinetic momentum. 



Notice that Eq.(3.17) is independent of t. In other words, the distribution is static. 



This is a natural result because we take the constant electric field which has existed 
from infinite past. Of course, each particles are accelerated by the electric field. The 
distribution, however, has already reached static state since the electric field has continued 
to create particles from infinite past. How the distribution reaches static state will be 
clarified in SecJ4.1[ 



Furthermore, one may notice that the distribution (3.17) depends on the longitudinal 



momentum p^. As mentioned in the introduction, longitudinal momentum dependence 
of the distribution conflicts with the boost-invariance of the electric field. In our calcula- 
tion, the boost-invariance is broken by the quantization procedure. To define a particle 
picture, canonical quantization is needed, which is done on an equal-time line (or at least 
a space-like line). That is to say, particle is global concept and its definition breaks the 
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boost invariance. To say more physically: before particles are created, the system with 
a constant electric field is boost-invariant along its direction. However, once a particle is 
created, it brings a special frame, namely center of mass frame [35] . 

Now the following question may arise: Is there a special preferred frame to describe 
pair creation? Suppose that observer A finds the distribution rip and observer B moves 
relatively to A along the longitudinal direction. A and B are on an equal footing, because 
the electric field is boost-invariant along the longitudinal direction and thus A and B feel 
the same electric field. Therefore, A and B must observe the same phenomena in each 
frame: A observes Up and B observes ripi, where p' is the momentum in the i?'s frame 
and each n(.) has the common functional form. However, if one supposes a phase space 



distribution is scalar with respect to Lorentz boost the distribution observed by B takes 
the different functional form from that by A: A observes rip, while B observes rip = n'^,. To 
say more specifically, if A observes that a particle is created with longitudinal momentum 
of (which is expected from the form of the distribution; Figjl]), then B observes this 
particle is created with non-zero momentum. In this situation, observer A looks like a 
special observer. 

Of course, it is not the case. Although the functional form of the distribution rip 
[Eq.(3.17)] is not boost-invariant along the longitudinal direction (i.e. n(.) ^ n'^.-j ), every 



observers observe the same distribution [modulo the ambiguity of particle picture in a 
background field]. Observer A finds the distribution rip based on the particle picture 
defined on an equal-time line of A's frame: rip = (0,in|a^ (t)ap(t)|0,in)^. On the 
other hand, observer B finds the distribution with the same functional form with A^s: 
ripi = (0, in'|ap,(t')ap'(t')|0, in')^^7-. This is because B uses the particle picture defined 
on an equal-time line of i?'s frame, which is different from ^'s. Thus, if one simply boosts 
the ^'s distribution to the i?'s frame, it is different from i?'s distribution because the 
boosted A^s equal-time line is different from the i?'s. That is to say, our distribution 



(3.5) is not scalar with respect to Lorentz boost unlike that in classical theory because 
quantization procedure requires specifying a Lorentz frame. If we take account of the 
fact that quantization procedure is not boost-invariant in the presence of a pair-creating 
background field, we can get the distribution, of which functional form itself is not boost- 
invariant, however, which is compatible with the boost-invariance of the background field. 



In Figjljthe pair distribution function (3.17) is plotted as a function of the longitudinal 



momentum p^ for several values of a = Hereafter we show all figures with dimension- 
less values which are scaled by y/eE [or y/eE^ in a non-steady field case, where Eq is an 
initial field strength]. 

The shape of the distributions can be understood roughly as follows: particles [an- 
tiparticles] are created with approximately 0-momentum, and accelerated by the electric 
field to get positive [negative] momentum. This process have been continued from infinite 
past and makes, roughly speaking, a step function-like shape of the distribution. 



The asymptotic forms of the distribution (3.17) at pz — > ±oo are 

Up e-2™ {p, +oo) 
{pz -oo). 



This is a matter of course by definition in classical theory. 
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Figure 1: Boson distributions 
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Figure 2: Fermion distributions 

The value at ~^ +oo agrees with which has been obtained in Secj2| where the asymp- 
totic WKB criterion was used. This agreement is consistent with the opinion made in 
Secjgl that we can see only particles with infinite longitudinal kinetic momentum using 
the asymptotic WKB criterion. 

One notable point seen in Figjl] is divergence of the distribution of massless particle 
at p = IpI = : 

, {r(!)rvi£ 



Hr- 



m=0 
p~0 



27rp 



(3.19) 



This divergence comes from the fact that the energy cost to create massless particles is 
0. If there is an electric field for any strength, it creates massless particles endlessly. It 
is adequate to treat a massless and 0-momentum mode as "field" rather than "particle" . 
Thus, if one interprets this mode as particles and counts their number, it is infinite. In 



this sense, the divergence (3.19) is of the same quality with that of the Bose-Einstein 



distribution of massless particle: 



n 



BE 



^LUp/T — ]^ m=0 p 



(3.20) 



p~0 



Therefore, if one watches only the behavior around 0-momentum, the distribution of mass- 
less particles created in the constant electric filed takes the same form as that of the Bose 
distribution. The electric field plays the role of a thermal bath, of which "temperature" 
is {r(|)}^v^/2vr. 
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Figure 3: Comparison between the distribution (3.17) and the "WKB interference term' 



(3.21), which is plotted with an arbitrary vertical scale factor. 



Other notable point seen in the distribution is oscillation with respect to pz- Its pe- 
riodicity is the order of \/~eE and slightly depends on a. This oscillation may be caused 
by interference between particles and antiparticles. Now, particles and antiparticles are 
created everywhere in whole space and accelerated by the electric field; particles are ac- 
celerated to the same direction with the electric field, while antiparticles to the opposite 
direction with the electric field. The WKB wave functions representing accelerated par- 
ticles and antiparticles are exp [i5'ci(t)] and exp [— 'i5'ci(t)] respectively, where Sdif) is a 



classical action of a particle under the electric field (2.20). Because the concepts of par- 
ticle and antiparticle are mixed by the Bogoliubov transformation in the electric field, 
particles and antiparticles can interfere with each other through the electric field. These 
interference may be represented as 



-iScl{t)\2 



2-2COS [2Sd{t)]. 



(3.21^ 



After changing the variable from t to pz through the classical equation of motion = eEt, 



we plot Eq.(3.21) together with the distribution (3.17) in Fig|3j The oscillatory behavior 
and its slight dependence on a are very well agreed, although a physical origin of Eq.(3.21 ), 
especially the reason why the relative phase between e^'^'^'*-*-' and e"*"^"'*^*^ is —1 is not clear. 
This agreement is limited in the region of > because particles are created with 
approximately 0-momentum and accelerated to get positive momentum. 

Finally, let us derive a source term of the Boltzmann-Vlasov equation with the as- 



sumption that the distribution function (3.17) satisfies the Boltzmann-Vlasov equation 

d 



d 

r.:np{x) + V ■ Vnp(x) + eE-—np{x) = S{x, p) 
at opz 



(3.22) 



Substituting the distribution ( |3.17[ ) into Eq.( |3.22[ ), we can get the source term as follows: 

- 2^2^— Re 



^(P) 




2{ujp-pz)- 



+ Y ^ [{ujp - Pzf + 2(cjp - pz)pz - m%] Re 



-ia-2 - 



(3.23) 
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(a) bosons (b) fermions 

Figure 4: The source terms 



where the arguments of parabohc cyhnder functions, (— e4*w^p^), are omitted. 



The 



source terms for bosons are shown in FigQa), in which oscillatory behavior appears. 
Similar oscillatory source term has been obtained also in Ref. [Ml 1381 l39t HO]. 

It will be convenient for a latter discussion to express the source term using the Bogoli- 

ubov coefficients. Differentiating (3;{t) = f^p4^p-eEt(^) - ^^n-eEt(^)l [Eq.(|3l2l) 



by 



dt 



and using the equation for j^™ [Eq.(2.5)], one can verify 
d 



dt 



OB* 83* 



dt 



dpz 



Pz 



ieE^t(3*(t) + eE^e'''^-'a^{t). 



2ul 



With this equation, we can derive the following expression for the source term: 



5(t,p) 



d 
dt 



2Re 



eE^Re [e-2-p*ap(t)/3p(t)] 



(3.24) 



(3.25) 



This is a general expression under a homogeneous and constant electric field. 



3.1.2 Collinear electric and magnetic field 

We derive a distribution function of created particles in a collinear electric and magnetic 
field in the same way as that in the pure electric field. 

There is a difficulty to define instantaneous positive and negative frequency solutions: 
The gauge potential at instant of time t = to, namely A{to,x.) = {—By,0, —Eto) is not 
pure gauge. It gives the constant magnetic field. Mode solutions under the gauge A{to, x) 
is not simple plane wave solutions: 

<!{to)(„\ — ^ ^.Tii^pt JPxX+i{p^-eEto)z I L I eB\^ 1 



±/^*«^(x) = — =6^*'^''*— e*P^^+*^P^-"^*o^^/— — ^Dn{Ti). (3.26) 



2ujp 271 \ 271 \ 71 



However, in spite of the presence of the magnetic field, we can use these solutions as 
instantaneous positive and negative frequency solutions. That is because the magnetic 
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field does not affect separation of positive and negative frequency solutions [e ^'^J'* and 
g+iojptj^ which means that a pure magnetic field cannot create particles. 



The use of the instantaneous positive and negative frequency solution (3.26) instead of 



(3.1) simply brings the replacement (2.48). Therefore, the resultant distribution function 



is (3.17) with the replacement a 



2eE ■ 



3.2 Fermions in a constant electromagnetic field 
3.2.1 Pure electric field 

In this subsection, we calculate a pair distribution function of fermions in the pure electric 
field using the instantaneous particle picture. 

The instantaneous positive and negative frequency solutions at time of t = to under 



the gauge potential = (0, 0, 0, —Et) is obtained from the free-field solutions (B.12) by 
the gauge transformation = —EtQ-. 



(3.27) 



The inner products between these solutions and the in-solutions ±^p'^g [Eq.(B.28)] are 



e 4 



ps 



qs' 



2VeE ^/ujp{ujp - p^) 

X 5s,/5^(p - q - eEto) 
ap{tQ)6ss'5^{v - q - eEto) 




Pz 




Pz 



X 6ss'S^{p - q - eEto) 
-/5;(toMss'5'(p-q-eEto 




Pz 




Pz 



(3.28) 



from which the relations 



.^-(X) = ap+eEto(^o)+<*:U.(^) -/5;+eEto(^o)-C+eE*o 



(to) 



(to) 



[X 



(3.29) 



are followecCB 



These expansions is definite by themselves and the additional condition (3.8 1 is not necessary, because 



the Dirac equation is a first-order differential equation and the associated inner product (2.60) does not 
contain time-differentiation. 
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The field expansions 



(3.30) 



and Eq.(3.29) yield the relations between the creation/annihilation operator of the in- 



particle picture and that of the instantaneous particle picture of t = to: 

«ps(to) = ttp(^o)ap-eEios + /^p(^o)&™p+eEtos ^ 



T,i). 



(3.31) 



Notice that this Bogoliubov transformation does not mix the two spin states and the 
Bogoliubov coefficients are common to both spin-| and | state because we use the same 
spin basis for both the in-solutions and the to-solutions. Thus, we can get the same 
distribution function from any of 

(0, in|aj,^(t)apt(t)|0, in) (0, in|a^^(t)ap|(t)|0, in) 

(0,in|6lp^(t)6_pt(t)|0,in) (0, in|6Lp;(t)6-pi(t)|0, in). 

Then, the pair distribution function at time t is 

--a 

e ■i 



AeEujp{ujp + pz 



m'^D^ia-i{-ei'' 




2 ^ ^ . / 2 




(3.32) 



This is independent of time with the same reason as that of boson. 



We plot the distribution (3^32) as a function of the longitudinal momentum pz for 

several values of a = ^ in Fig|2 

The fermion distribution shows little oscillation compared with boson. This is because 
of Pauli's exclusion principal, which makes collective behavior of fermions "particle-like", 
while bosons behave like as "field-like" . This "particle-like" nature of fermions enhances 
the step function-like forms of the distributions. Especially in the case of a = 0, mass- 
less fermions are created as many as possible under the restriction of Pauli's exclusion 
principal, therefore, the distribution becomes the step function exactly: 



_ \Pz \ +Pz 
a=0 2\pz\ 



(3.33) 



On the other hand, bosons which have "field-like" nature interfere during acceleration 
and their distribution shows oscillation. 

Asymptotic behavior of the distribution (3.32) at p^ —>■ ±oo are just the same with 
those of bosons: 



-2na 



rir 







{p 
{p 



-oo ) 
-oo) 



(3.34) 
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This is because, as already noted on Sec, 2.2, high energy hmit {pz — > ±00) bhnds spin 
statistics. Meantime, low momentum region greatly reflects spin statistics. 



The source term obtained by substituting the distribution (3.32) into the Boltzmann- 



Vlasov equation (|3.22|) is 



S{p) = eE 



UJf, 



+ eE 



e 2' 




Pz 




(3.35) 

which is plotted in FigQb). The source terms for fermions have a peak at Pz = 0. 
Especially if a = 0, the source term becomes 



Sip) = eE5{pz). 



a=0 



(3.36) 



This delta function source term has been introduced phenomenologically in Ref . [351 136j*^^[ 
On the other hand, when a 7^ 0, the peak at = is broadened by quantum fluctuation. 

Finally let us derive a general expression of the source term under a homogeneous 
electric field. Using expressions 



[mFi - {p, + tpy)r2 + {±00, - p,)r3] sPi°\x), 



and 



+ij''^{x) = Afp [mFi - (p, + zpy)T2 + (ido -Pz- eEt)T,] +0j(t) 



(3.37) 



(3.38) 



in which +(/'p(t) is the positive frequency solution of "Klein-Gordon equation" (B.6) and 
Afp is a normahzation constant, one can get the equations 



m - 



^Jujp{u}p + p^) 



.^p+Pz)Mp-eEt-P. 



p_eEt - Pz+<Pp-eEt) + "^T+0p-eEt) 



zWp-eEt) ~ ""^T+^p-eEi) 



(3.39) 
(3.40) 



Differentiating (3.40) with respect to t and using Eq.(B.6) yield 



Z J. QJf 



(3.41) 



Then, the source term 



5(t,p) = -|/?p(t)r = eE^Re[e- 



a, 



(3.42) 



is derived. This will be useful for a latter discussion on a polarization current. 



*12 



However, the argument of delta function is not momentum but rapidity in these reference. 
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3.2.2 CoUinear electric and magnetic field 



Similarly to the boson case, we use solutions under the gauge potential A(to,x) = 
{—By, 0, —Eto) as instantaneous positive and negative frequency solutions at time of 
t = to. They are given as follows: 



X 



ps 



x e 



-iEtoz 



--hi), 



(3.43) 



in which -t^/^r, is solutions under the gauge A = {—By, 0,0) [Eq.(B.31)]. The pair dis- 



tribution function can be obtained by taking the inner product between the in-solutions 



(B.23) and the instantaneous solutions (3.43). In the same way as Sec, 2.2.2, the longitu- 



dinal magnetic field results in the replacement (2.70). 



4 Non-steady electric fields 

In the previous section, we revealed the momentum distribution of created particle using 
the instantaneous particle picture. However, if one integrates the distribution (for example 



Eq.(3.17)) with respect to the momentum, it diverges although it represents the total 



particle number per unit volume. This is not a defect of the instantaneous particle picture 
but a natural result of the constant electric field; particles have been continued to be 
created from infinite past and their number has already reached infinity. 

To get a finite particle number per unit volume, we must treat a non-steady case in 
which an electric field strength is up to some time. It will give us the information how 
the distribution evolve in time. 



4.1 Sudden switch-on 

As the simplest example of non-steady electric fields, we first deal with the "sudden 
switch-on" electric field: 

E ={0,0, Ee{t- to)). (4.1) 

Although realizing this space-time configuration is impossible due to causality, we use it 
because it is the simplest model revealing the dynamical process of pair creation and a 
distribution function in this electric field can be calculated analytically. 



The gauge potential giving the electric field (4.1) is 



-Et {t > to) 
-Eto {t<to). 



(4.2) 



We show only the derivation of the boson distribution function since the fermion 
distribution can be obtained in a similar way. 

In this case, we do not need to use the asymptotic WKB criterion because the electric 
field does not exist before the time to. The in-solutions at t < to is obvious: 



(to) I 



(for t <to 



(4.3) 



where ±fp°^ is given by Eq.(3.1). One can get the in-solution at t > to by solving the 
field equation with Eq.(4.3) as an initial condition at t = to- However, we do not need to 
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re-s olve t he field equation. We c an co nstruct the in-solutions at t > to using the result in 

Sec, 



3.1.1 



follows: 



Rewriting in Sec 3.1.1 as ±f^, we can obtain the in-solutions at t > to 



+/^(a;) = a;(to)+/J_eE*o(^) - /^;(^o)-/jr_eE*o(^) 



(for t > to 



(4.4) 



where ap(to) and /?p(to) are defined by Eq.(3.12). One can confirm that this solutions 
satisfy the initial condition at t = to by substituting Eq.(3.13) (which holds only at instant 



of time to) into ±fp. We can get the relation between the in-solutions (4.4) and the 
instantaneous positive and negative frequency solutions at time ti(> to) by substituting 



Eq.(3.13) with the replacement to — > ti into ±fp of Eq.(4.4) 



,/j^(x) = Ap(to,ti)+/, 



(*i) 

p+eEAi 



B;{to,h)jl%l^^,ix) 



where 



^p(i^O,^l) = «p(to)«p+eEAi(tl) - /3*(to)/9, 



'p+eEAt(tl 



(4.5) 



(4.6) 



Bp{to,ti) = Q;p(to)/3p+eEAt(tl) — /3p(to)a;p+eEAt(tl) 

and At = ti — to. 

To define an in-particle picture and the associated vacuum, we expand the field oper- 



ator with which is given by Eq.(4.3) and (4.4): 



m = / dpy'^{x)a^ + J^ix)b 



intl 



(4.7) 



The Bogoliubov relation between the creation/annihilation operator of an in-particle and 



that of a ti-particle can be found by substituting Eq.(4.5) into Eq.(4.7): 



0(x) = / d'p [ap(ti)+/(*^)(x) + &i(ti)_/^^)(x)] 



ap(ti 



A 
A* 



p-eEAt(to, ^l)'^p-eEAt + -Bp-eEAt(to, ^l)^-p+eEAf 



p-eEAt 



(to,ti)6 

-p+eEAt "T -^p-eEAt 



+ B* 



(to,ti)<- 



(for ti > to). 



(4i 



(4.9) 



p-eEA< 



Now we can calculate the pair distribution function at time ti(> to) as follows: 



(27r) 



np(At) = (0,in|aT,(ti)ap(ti)|0,in) ^ - \^,x..|„_p 

= \Bp-eEAt(tQ,ti)\'^ 



(0,in|6lp(ti)6^p(ti)|0,in) 



(2vr) 
V 



(4.10) 



This depends on time through only At, because ap(t) and Pp(t) depend on t through only 
the phase factor e^*'^*'*. Of course, np(At) = when ti < to. 

The distribution of fermions can be derived in a similar way. The result is the same 



instead of Eq.(3.12) 



form as that of bosons [Eq.(4.10)], however ap(to) and (3p{to) have been given by Eq.(3.28) 
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(c) fermion a = 0.3 (d) fcrmion a = 0.01 

Figure 5: The time evolution of the longitudinal momentum distributions 

The time evolution of the longitudinal momentum distribution, n.p{t) with fixed a, and 
that of the transverse momentum distribution, Tip(t) with fixed p^, are shown for both 
bosons and fermions in FigjS] and Figj6j 

The longitudinal momentum distributions rise up first around 0-momentum after the 
field is switched on. As already expected in SecjSj only particles with approximately 
0-momentum can be created. However, not strictly but broadened around due to 
quantum fluctuation. 

Comparing the two longitudinal momentum distributions for different values of a, we 
can notice that the time when the distribution peaks first is earlier for smaller a. This 
is to say, the lighter the mass or the stronger the electric field, then the shorter time the 
electric field takes to create particles, or for fixed eE and m, particles with pt = are 
created first and those with higher px are later. 

After particles are created with about 0-momentum, they are accelerated by the elec- 
tric field according to the classical equation of motion: pz = eEt + const. This process 
makes triangle shaped distributions in t-p^ space. On the other hand, the transverse 
momentum distributions take a saddle-like shape in t-p^ space. This is because the trans- 
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(a) bosons 0.3, pjV^^ = 1.45 (b) fermions 0.3, pjV^^ = 0.77 

Figure 6: The time evolution of the transverse momentum distributions 



verse momentum distribution can be approximated as a Gaussian at any time and settles 
in a fixed distribution after sufficiently long time; especially e~^™ in the limit of pz oo. 

Both of the longitudinal and the transverse momentum distribution of bosons show 
oscillations. These oscillations are the same as that has seen in the distribution in the 
constant electric field [Figjl] 



One can recover the distribution in the constant electric field [Eq.(3.17) or (3.32)] 
taking the limit of At +oo in Eq.(4.10). Therefore, the use of the asymptotic WKB 



criterion in a constant electric field is consistent with the treatment in this section, where 
the in- vacuum is defined in absence of an electric field and the asymptotic WKB criterion 
is not used. 



Let us comment on a symmetry of the distribution (4.10). The distribution np{t) 



is symmetric by reflection with respect to the axis pz = i^^Et in pz space; 'T-ieEt+pl^) 



ni. 



Et-p 



(t) for any t. This symmetry originates in the time reversal symmetry: (i) Because 



the argument (t) of a^if) and f3p{t) is irrelevant to np{t), 



-t) 



(4.11) 



holds formally for t > 0. We put^'on n_p(— t) because originally Eq.(4.10) gives np{t) only 
for t > 0. (ii) Seeing the derivation of Eq.(4.10), one can notice that ^p(t) of Eq.(4.10) 
for negative t represents the distribution in the time-reversed system, namely the system 
with the boundary condition ^p(O) = in the electric field which is switched on at infinite 
past and switched off at t = 0. Thus, 



?2_ 



-t) = Tn_p(-t) 



(4.12) 



is satisfied, where T is a time-reversal operator, (iii) Because this system is invariant 
under time reversal, 

Tn_p(-t) = np{t) (4.13) 



holds. Eqs.(4.11 ),(4.12) and (4.13) yield the symmetry in question: np{t) = rieEt-pit). 



Finite time imposition of the electric field offers us finite expectation values of the 
total particle number density, the charge current density and so on. In Fig|7|^a) we show 
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(c) the polarization current of bosons 




■ « = 0.3 
a = 0.01 



4 6 

(d) the polarization current of fermions 



Figure 7: The time evolution of the total particle number density and the current density 



for several values of a 
states, e = 1. 



2eE 



. For fermions, these values have been averaged over the spin 



the time evolution of the total particle number density N = J yf^^rs^p for several values 



of a = In Fig,7|^b),(c),(d) the current density, which is separated into a conduction 
current and a polarization current is plotted. For fermions, the values averaged over the 
two spin states are shown. 



The Noether current associated with U{1) symmetry in Lagrangian (2.1) is 

= ie [(f)\d^ + teA^')(f) - {d" - ieA^')(j)'^ ■ 0] . (4.14) 
Taking the in-vacuum expectation of the z-component of this current operator with the 



ti-field expansion (4.8) and the gauge potential at instant of time t = ti. A = —Eti, we 



get the mean current density 
j,{h) = (0,in|i|0,in) 



e / (fp / d^q 



Pz + Qz 



:(0,in| [, 



al{ti)a^{ti) 



(27r)3cjp"^P' J (27r)3ajp 

f d?p Pz /, N , o /" Q(+ \ 



bosons. 



(4.15) 
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Figure 8: Schematic picture of a tunneling process 



The first term is a conduction current, which is caused by movement of real particles and 
the second is a polarization current, which originates in pair creation process. 

A schematic picture of a tunneling process, FigjS} may help us to understand a meaning 
of the polarization current. Figure [8] itself is applicable to only fermions, however, the 
following explanation is also relevant to bosons. Suppose that a pair of a particle with 
momentum p and an antiparticle with momentum — p is created. The distance between 
them is 2^ when they become on-shell. Thus, this process can be understood as a 
creation of an electric dipole 2e^. The polarization current is induced by the variation 
of an electric dipole and now the creation rate of a dipole is ^np(t) = S{t, p). Therefore, 
the variation rate of the total electric dipole, or the polarization current is 

^r'W=/(^2e^S(«.p). (4.16) 



which is identical to the second term of Eq.(4.15). 



For fermions, the z-component of the mean current density is 
jz{t) = -(0,in|eV'7^?/'|0,in) 



(4.17) 

fermions , 



of which last expression is identical to that of bosons except spin summations. 

The polarization current has a pulse just after the switch-on. This is because the 
particle number undergoes a rapid increase from the initial state with no particle. Some 
time after the switch-on, the increase of the particle number gets equilibrated to a fixed 
rate. Therefore, the polarization current settles in a constant value and the conduction 
current becomes proportional to time. 

Although we have the analytic form of the distribution functions, carrying out the 
integral analytically in the total particle number density N = J 7^^^p(^) and the current 



density (4.15) is impossible and the parameter dependence of these values is not clear. 



Thus, we model the distribution function in a simple form 

np{t) = e-2™e(p,)^(eEt - p,), (4.18) 
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Figure 9: The total particle number density and the current density of fermions in the 

2 

collinear electromagnetic field (a = ^ = 0.01, e = 1). The values averaged over the spin 
states are shown. 



and investigate the parameter dependence. This approximated distribution is exact for 
fermions with a = and valid for other cases except the time just after the switch-on and 
bosons with a small a. Eq.(4.18) corresponds to the delta function source term: 



S(t, p) = eEe-'^^d^p, = O)^(t) 



(4.19) 



With the help of Eq.(4.18) and approximations which are valid under m <C E, the 



integrals in the total particle number density and the current density can be easily done: 



N{t) ~ ; ; e-^t e{t) 



(2vr) 
2e{eEf 

e 



(27r) 



t 9{t) 



, , e(eE)^/^ ™2 



(4.20) 
(4.21) 
(4.22) 



These are common to bosons, fermions with spin-f and with spin-| mode. 

The approximated values ( |4.20 ), (4.21) and (4.22) succeed in reproducing the slope 
of the linear increase of N{t) and jz°'^^{t) obtained by the numerical calculation and the 
constant value that j^°^(t) approaches, except (i) the time just after the switch-on, in 
which the system remembers the information of the initial state, on the other hand, after 
sufficiently long time from the switch-on, the pair creation process settles in a steady 
state, which is well expressed by the approximation (4.19), and (ii) bosons with a small 
a because its distribution has large peaks on the Pz = and pz = eEt lines which is not 



expressed in the approximation (4.18). These approximated expressions will be useful to 



analyze a time scale of back reaction (SecjS 

In Fig|9]the time evolution of the total particle number density and the current density 
of fermions in the collinear electric and magnetic field are shown. (We do not plot those 
of bosons because they are strongly suppressed by an intense magnetic field.) As already 



noted in Sec. 2.2.2 the longitudinal magnetic field increases balk quantities of fermions. 
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In the collinear electromagnetic field, the approximated particle number and currents 

are 



3 



Nit) 

cond 



2e'^EB 



it) 



(27r)2 
2e^EB 

N3/2 



e 



1 - e 



7rm 

e 



-t e{t) 



t d{t) 



(4.23) 
(4.24) 
(4.25) 



n=0 



where 



1 for bosons 

for fermions with spin- | (4.26) 

2 for fermions with spin- J, . 

Also in this case, the approximated values quite well agree with the exact ones with the 
same exception as before. Eq.(4.23) tells us that the longitudinal magnetic field increases 
the particle number of fermions with spin-| approximately in proportion to B. 

So far, we have employed the in-vacuum as an initial state, in which case there is no 
particle before the electric field is turned on. We can also use the state where particle 
pairs are condensed, in which case there are initial particles and they evolve under the 
influence of the electric field after the switch-on. 

The state with an arbitrary initial distribution /(p) is constructed as follows: 



|/)=nA/'p^/'exp ^F(p)a-t6L^t 



|0,in). 



(4.27) 



where -F(p) and the normalization constant A/'p are related with the initial distribution 
/(P) as 



I^(P) 



/(P) 



+ : bosons 
— : fermions 



(4.28) 



1±/(P)- W-fli/C)}" 

One can confirm that the expectation of the in-particle number operator with the state 
I/) is surely equal to the initial distribution /(p): 



(27r 



(2vr) 
V 



/(P)- 



(4.29) 



The distribution at time t with the state |/) is 



(2vr) 



np(t) = (/|at(t)ap(t)|/) ^ 



(/l&lpW&-pWI/)^ 



[1 ± 2np(At)] /(p - eEAt) + np(At) 

[1 ± 2/(p - eEAt)] r2p(At) + /(p - eEAt), 



(4.30) 



in which rip (At) is the distribution with the state |0, in) [Eq.(4.10)]. This distribution can 
be interpreted as the summation of the initial distribution accelerated according to the 
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(a) bosons, a ~ 0.01 (b) fcrmions, a = 0.01 

Figure 10: The longitudinal momentum distributions with the initial condition /(p) 



classical equation of motion /(p — eEAt) and the usual distribution of created particles 
rip (At), however, there is a quantum correction; Bose enhancement or Pauli blocking 
factor 1 ± 2/(p — eEAt) [or one can consider 1 ± 2np(At) are that factors.] 

The Bose enhancement factor tells us that if there are initial particles, the pair creation 
at the point where the initial particles are located in momentum space is enhanced. 
On the other hand, for fermion system, initial particles suppress the subsequent pair 
creation. Especially, if np(At) > 0.5, an initial distribution reduces particle number; 
np(At) < rip (At). As an illustration of these phenomena, we show the time evolution o: 
the distribution (a = 0.01) with an initial condition of /(p) = O.Se"^^/^^"*"^^^ in Fig, 10 



Actually, the boson distribution increases more than double when the initial particles 
merge into the region of rip (At). Meanwhile the fermion distribution gets dented because 
?7,p(At) > 0.5 holds for a = 0.01. These effects will be conspicuous in the evolution of a 
distribution with back reaction (Sec jsj). 

If one sets |F(p)p = e"'^''/^, Eq.(4.30) describes the Schwinger mechanism at finite 
temperature. Our result is slightly different from those in Ref . |l9l [56] , where the Schwinger 
mechanism at finite temperature is studied by means of the density operator. In Ref. |l9l 
[56] , thermal distribution is set to be free from acceleration by the electric field. 



4.2 Damping electric field 

In the previous subsection, we have dealt with the electric field which is kept constant 
after the switch-on. However, in most of realistic situation, an electric field is damped 
away after sufficiently long time, because (i) its energy are taken away by particles due 
to pair creation and subsequent acceleration, and (ii) energy of particles would dissipate 
owing to collision among particles. The reason (i) will be treated in the next section, 
however, (ii) is not considered in this article. Furthermore, in heavy-ion collisions, there 
is another mechanism which makes an electric field decrease; classical time evolution of 
gauge fields following expansion of the system [521 [66] • 

Instead of studying these effects directly, we treat a damping electric field, which is 
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(a) bosons a = 0.3 



(b) bosons a = 0.01 





(c) fcrmions a = 0.3 



(d) fermions a — 0.01 



Figure 11: The time evolution of the longitudinal distributions in the damping electric 



field (4.31) with r = l.O/v^- 



controlled by hands. As an example, we take an exponentially decreasing electric filed: 

(4.31) 



^ ' ^ Eoe^'l- > 0) ^ ' 



(t<0) 
EoTe'^/^ (t > 0). 



We have solved the Klein-Gordon equation and the Dirac equation numerically. 

A generalized expression of the Bogoliubov coefficients under the gauge A{t) = (0, 0, A^{t)) 

is 



ap{t) 



(4.32) 



for bosons and 

ap(t) = , 

13* it) = e-'"P* _ 



{Ujp - P.)(n0;,\,A - P.+<+eA) + ^T+0;.\eA) 



(4.33) 
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(c) the polarization current 
(b) the conduction current 

Figure 12: The time evolution of the total particle number density and the current density 



in the damping electric field (4.31) with r = ^/cEq^c = 1. For fermions, these values 



are plotted after averaged over the spin states. 



for fermions, where i(/>p (t) are the solutions of Eq.(A.l ) for bosons or Eq.(B.6) for fermions 
with the initial condition i(/>p (t) = ^1 — ^T^i^pt f^j. t < respectively. Using these equa- 



tions, we can calculate the distribution function, the total particle density and the charge 
current. 



The results of the numerical calculations are shown in Figs, 11 and 12 The distribution 
functions evolve in time naturally as one can expect from the results in the previous 
subsection (Figjsj): (i) rise up in a similar way to Figjs] (two peaks structure in FigjTTj^a) 
and (b) are also seen in Figjsj^a) and (b) just after the switch-on), and (ii) accelerated to 
get higher momentum, however, (iii) pair creation and acceleration is stopped after the 
electric field is damped away and there left a steady distribution. 

Before the switch-on and after the electric field is damped away, our particle picture is 
absolutely reliable owing to absence of the electric field. Furthermore, the instantaneous 
particle picture gives us the proper description of the time evolution of the system during 
intermediate time even in the presence of the electric field. During the evolution, there is 
no pathologic behavior in the distribution or the total particle density unlike the case in 
cosmological particle creation [T3] . This observation as well as the results in the previous 
subsection would support the validity to use the particle picture associated with the in- 
stantaneous positive and negative frequency, and to study time evolution of a distribution 
in the presence of an electric field. 
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(c) fermion a = 0.3 (d) fermion a = 0.01 



Figure 13: The time evolution of the longitudinal momentum distributions with back 
reaction (e = 1). 

5 Back reaction 

So far, we have controlled the electromagnetic fields by hand and neglected back reaction; 
an electromagnetic field which is generated by particles and antiparticles created by an 
original electric field. To take account of back reaction, we treat an electromagnetic 
field as a dynamical variable. However, we does not quantize it and treat it as a classical 
background. This approximation may be valid for a strong electromagnetic field, quantum 
fiuctuation of which is negligible. 

Now we use the Lagrangian density 

C = ~\f^uF^'' + {d^ - zeA^)0t . (9^ + zeA^)0 - m^^, (5.1) 

in which F^^ = d^Ay — d^A^. The equations of motion followed from this Lagrangian are 
the Klein-Gordon equation 

[(9^ + leA^f + m^] (P{x) = (5.2) 

and the Maxwell equations 

V-E=jo, ^-VxB = -j, (5.3) 
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Bosons, a = 0.01 
Fermions, ci = 0.01 

Bosons, a = 0.3 
Fermions, a = 0.3 
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(a) the total particle number density 




(b) the current density (c) the electric field 

Figure 14: The time evolution of the total particle number density, the current density 
and the electric field (e = 1). For fermions, the values averaged over the spin states are 
shown in (a) and (b). 



in which is given by Eq.(4.14). The current operator is replaced by its expectation 
value (0, in|j^|0, in) because we treat the gauge field as a classical field. A homogeneous 
initial electric field directing along the z-axis reduces the Maxwell equations into the single 
equation 



dt 



dt^ 



-3z 



(5.4) 



where the 2;-component of the mean current jz is given by Eq.(4.15) for bosons or Eq.(4.17) 
for fermions. 

We have solved the Klein-Gordon equation (5.2) [or the Dirac equation] and the 
Maxwell equation (5.4) numerically with initial condition of the "sudden-switch-on" : A'^ = 
0, = for t < and A^ = 0,A^ = —Eq at t = 0. We calculate the distribution function 
and the charge current with the help of Eq.(4.32) or (4.33). 



The results of the numerical calculations are shown in FigsjT4]and[T3] The longitudinal 
momentum distributions oscillate in momentum space and also the current and the electric 
field show oscillation. Similar oscillatory behavior have been obtained in earlier works 

ESI El EZl SD 112], although the employed methods are distinct. 

These oscillations can be understood as usual plasma oscillation. First, particles [an- 
tiparticles] are created with approximately longitudinal momentum. They are accel- 
erated to get positive [negative] momentum, and generate positive conduction current. 
Then, the electric field is reduced by the current following the Maxwell equation (5.4). 
After a while its strength reaches 0, at which time acceleration of particles is stopped, 
however, the current is still positive because particles have already gotten positive mo- 
mentum. Thus, the electric field continues to decrease through to negative value (the 
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direction of the electric field is flipped). Then, particles start to be decelerated and finally 
their momentum become 0, at which time decrease of the electric field is stopped. How- 
ever, particles continue to be decelerated and get negative momentum, which generate 
negative conduction current, because the direction of the electric field is negative at that 
time. Then, the negative conduction current increases the electric field and the direction 
of the electric field is flipped again. These processes are repeated and result in plasma 
oscillation. 

The oscillation is faster for (i) a = 0.01 than a = 0.3, and (ii) fermions than bosons. 
These are because the current is larger for (i) smaller value of a, and (ii) fermions than 
bosons due to spin degrees of freedom. 

Although plasma oscillation itself is a purely classical dynamics, also quantum effect, 
namely Bose enhancement or Pauli blocking plays an important role during the time 



evolution. As explained in Sec, 4.1, particles of boson [fermion] enhance [suppress] the 
pair creation at the point where they locate in momentum space. These effects are more 
conspicuous for smaller a, because the distribution depends on a roughly as e 



-2-Ka 



see 



Eq.(4.18)] and the effects of Bose enhancement and Pauli blocking are stronger for larger 



value of the distribution [see Eq.(4.30)]. 

When the distribution crosses the line of = in momentum space, the distribution 
of bosons with a = 0.01 shows dramatic increases. Because pair creation can occur only in 
the vicinity of the Pz = line, pair creation is enhanced when old particles come back on 
the line of = in momentum space. Furthermore, the electric field takes its maximum 
value when the distribution crosses the Pz = line, because the current density is at 
that time. This also enhances pair creation. 
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(b) the current density (c) the electric field 

Figure 15: The time evolution of the total particle number density, the current density 
and the electric field for fermions with a = 0.01, e = 1. The values averaged over the spin 
states are shown in (a) and (b). 
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On the other hand, the distribution of fermions with a = 0.01 shows quit different 
behavior from that of bosons. It repeats sharp increases and decreases when it crosses 
the hne of pz = 0. This is because, as seen in Sec, 4.1, initial particles with Hp > 0.5 cause 
pair annihilation. 

Furthermore, we can apply a magnetic field parallel to the electric field. In this 
homogeneous system, a magnetic field is unchanged by back reaction. Thus, the effect 



of the magnetic field can be taken into account by the replacement (2.48) for bosons or 



(2.70) for fermions. We show the time evolution of the total particle number density. 



the current density and the electric field for fermions in the longitudinal magnetic field 



in Fig, 15 In the magnetic field, also the total particle number density shows oscillation. 
This is because of Pauli blocking and suppression of non-zero Landau level: (i) the fermion 
distribution with small a repeats sharp increases and decreases [see Figjl3|(d)], and (ii) the 
lowest Landau level, that is the smallest a mode mainly contributes to the total particle 
number in a strong magnetic field, while in the absence of a magnetic filed, larger a modes 



e.g. Fig, 13 ^c)] also contribute. 



The magnetic field enhances the current density of fermions by the mechanism ex- 
plained in Sec, 2.2.2 With back reaction, enhancement of the current makes the time 
evolution faster. This mechanism may be relevant in the time evolution of the Glasma 
flux tube. 



Energy balance between particles and the electric field is shown in Fig, 16 for fermions 
with d = 0.01, e = 1. Total energy conservation is automatically guaranteed in our particle 
picture and associated regularization method [3l] (see also [36j). The expectation value 
of the energy density operator, which is regularized by normal ordering is 2 J ■^:^ujpnp{t). 
Then, the time-derivative of total energy e is 

-ujpnp{t) 



de 
~dt 



d 
dt 



-E^ + 2 



E^ + 2eE 
dt 



e(^+ + jP°' 



i2n) 

d^P Pz 
(2vr) 



np{t) +2 j 



d^P . 
^-p5(t,p) 



(5.5) 



V dt 



in which 



dt 



dt ^ dp. 



and the Maxwell equation (5.4) are used. 
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Figure 16: Energy balance 
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Before closing this section, let us estimate a time scale of back reaction. As a time 
scale, we define the time when the electric field first reduces its strength to as ty. 
Although we have obtained ty by the numerical calculation, we estimate it analytically to 
see its dependence on parameter. As zeroth-order approximation, we employ the modeled 



expression of the conduction current without back reaction (4.21): 



J 



cond 



7rm 
eE. 



(5.6) 



in which the electric field is fixed to the initial value and Nd is introduced as a number of 
inner degrees of freedom: spin, flavor and color. Neglecting the polarization current, we 
substitute Eq.(5.6) into the Maxwell equation (5.4) and solve it. Of course, an oscillatory 



behavior cannot be obtained by this approximation, however, ty can: 



tv 



(27r) 







(5.7) 



Eq.(5.7) is an underestimate because we use the current without back reaction. Actually, 
Eq.(5.7) gives a smaller value of ty than the result of numerical calculations. However, 

J I 2 

their discrepancy is less than factor 2 and Eq.(5.7) correctly describes the a = and 
boson-fermion dependence of ty obtained by numerical calculations. Thus, we can safely 



use Eq.(5.7) for order estimation of t\ 



From Eq.(5.7), we can notice that the stronger the initial electric field, the shorter time 



it takes to vanish. This is because the conduction current (5.6) increases quadratically 
with increasing Eq. 

In the presence of the magnetic field parallel to the electric field, ty for fermions is 
estimated as 



tv 



Bo 



which is a decreasing function of the magnetic field strength Bq. 



(27r) 



tanh7r|a 



(5.^ 



6 Summary and discussion 

Our aim has been to make the physics in the Schwinger mechanism clearer. In particular, 
we have explored the momentum distributions of created particles and its time evolution. 
It has been accomplished by defining the particle picture at each time using the instanta- 
neous positive and negative frequency solutions. In the homogeneous and constant electric 
field, the static distributions have been obtained. We have shown that the longitudinal 
momentum dependence of the distribution arises from the boost-invariant electric field 
through the quantization procedure. Although this situation, namely the homogeneous 
and constant electric field is artificial, our consideration on the boost-invariance may have 
practical importance because electric fields formed just after heavy-ion collisions are also 
(approximately) boost-invariant along the longitudinal direction. 

The dynamical view of pair creation has been acquired by dealing with the non-steady 
electric fields. After the electric field is switched on, the distribution rises up around 
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777.^ I 

0-momentum. The time the distribution peaks first is earher for smaller a = ^^^'^ . 
That is, for fixed m and eE, particles with pt = are created first and those with 
larger px are later. In particular, the time to create a particle with a = is 0. After 
sufficiently long time, the transverse momentum distribution settles in a Gaussian-like 
form (Fig|6]). On the other hand, in the longitudinal direction, particles are accelerated 
by the electric field according to the classical equation of motion pz = eEt + const (Figj5]). 
Also the time evolution of the total particle number density and the current density have 
been computed. They are suppressed by the longitudinal magnetic field for bosons and 
enhanced for fermions. 

Taking back reaction into account, we have shown how the coherent initial electric 
field decays into particle pairs. The electric field loses energy owing to creating and 
accelerating particle pairs. We have estimated the time ty when the electric field first 
reduces its strength to 0. In our calculation, however, the electric field does not vanishes 
away after ty- It continues to exchange energy with particles and shows plasma oscillation. 

This plasma oscillation would be disturbed in reality by collision among particles, 
effects of which are not included in our calculation. We have treated the gauge field as 
a classical field. In this approximation, only particles and antiparticles having approxi- 
mately 0-momentum can be created or annihilated by the classical electric field. However, 
in reality, particles and antiparticles having non-zero momentum can collide and exchange 
their momentum through a quantum gauge field. These effects are expected to become 
relevant at t > ty? because the time evolution at early time t < ty is mainly controlled by 
the strong classical electric field. Perturbation may serve to treat these effects. However, 
usual perturbation based on asymptotic fields is no use. Real time formalism is needed. 
It is hoped to investigate this problem. 

To check whether this decay mechanism of the initial electric field due to back reaction 
is relevant to heavy- ion collisions or not, let us estimate ty with parameter motivated by 
the Glasma; m = and cEq ~ ~ 1 GeV^, where Qs is a saturation scale. Then, 



Eq.(5.7) yields 

tv 



If one simply assumes that quarks and gluons in a non-Abelian electromagnetic field can 
be treated as fermions (spinor particles) and bosons (scalar particles), respectively in an 



Abelian electromagnetic fielc , then Nd = 2 x Nj x + 2 x [N^ - 1 - {Nc - 1)] 



2x3x3-1-2x6 = 30, where iVg — 1 is the dimension of the Cartan subspace of SU(Ai"c). 
Thus, if > ^^(1); tv can be the order or less than 1 fm. This result would verify 
the importance of the Schwinger mechanism at the initial stage of heavy-ion collisions, 
although our estimation is rather rough. Of course, to discuss thermalization of the 
system, we must take a quantum gauge field into account. 

Furthermore, we have shown that the longitudinal magnetic field speeds up the time 
evolution of the fermion system and makes ty smaller. In the case of cEq ~ cBq ~ Q| ~ 



1 GeV^, ty gets smaller than Eq.(6.1) by the factor ~ y/n. 

In this paper, we have dealt with the spatially homogeneous electromagnetic fields. 
However, the fiux tube has a finite size in both the transverse and the longitudinal direc- 



*i3pQ]. quarks this assumption is correct [5], however, for gluons not correct because of intrinsic magnetic 
moment |67j . 
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tion, and expands to the longitudinal direction. Thus, for application to particle produc- 
tion in heavy-ion collisions, we need to consider finite size effects. Effects of transverse 
finite size have been studied in Ref. [271 1281 129] . however, the flux tube has been treated 
as a static object. In reality, the flux tube evolves dynamically by back reaction. In 
particular, back reaction to a magnetic fleld can be taken only if one treats a non-uniform 
fleld. Effects of longitudinal boundaries, especially in the case that boundaries move with 
the speed of light is also open to argument. Furthermore, consideration for interaction 
and correlation between flux tubes may have importance. Our present study would give 
the basis for further investigation of these effects. 

Note added in proof: After the completion of this work, the auther become aware that 
the time scale of back reaction has been studied by S.P.Gavrilov and D.M.Gitman [69] as a 
consistency restriction on maximal electric fleld strength, and the condition corresponding 
to Eq.(5.7) of us has been obtained. 
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Appendix A First-order Klein— Gordon equation 



In this appendix we briefly review a method to convert the Klein-Gordon equation 



\2 I 21 

+ m J 







into a flrst-order differential equation with respect to time-derivative |68j . 
We introduce the 2-component expression 



+ ^(9o + 2eAo) 



Then, the Klein-Gordon equation (|A.1|) can be rewritten as 
(9o + ieAo)$ 



2m 



(V — ieA) {ia2 + cr-^) — ima^ 



(A.l) 



(A.2) 



(A.3) 



which is flrst order with respect to time-derivative (cTj are the Pauli matrices). The 2- 
component expression reflects the degrees of freedom of particle and antiparticle. 
The inner product 



d'^x 



(pi ■ {do + ieAo)(f)2 - {do - ieAo)(f)\ ■ 02 



is expressed by <I>(x) as 



2m / d^x^{a3^2, 



(A.4) 
(A.5) 



which does not contain time-derivative explicitly unlike Eq.(A.4). 
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Appendix B Second-order Dirac equation 

We explain a method to convert the Dirac equation 



{B.l] 



into a 2nd order differential equation and solve it in the presence of a background elec- 
tromagnetic field [H]. In the following, we employ the 7-matrices in the standard Dirac 
representation: 



7 



1 
-1 



a' 
-a' 



Inserting the equation 



^(x) = [-f^iid^ - eA^) + m] Z[x) 



into (B.l) yields the equation for Z{x): 



le 



Z(x) 



0, 



(B.2) 



(B.3) 



(B.4) 



where Z{x) is a 4-component spinor and F^^, is the field strength: F^j, = d^A^, — d^A^. 
The term ^'^^'^^F^^ represents spin-electromagnetic field interaction. 

Because ^'j^'j'^Ffj^^, has 4 x 4-matrix structure with respect to spinor indices, it has 
four independent eigenvectors Ti{i = 1,2,3,4). We can expand Z{x) by Fj as 



Z{x) = J2<P\^)^i 



(B.5) 



Then, each obey the Klein-Gordon- type equation 

[{d^ + leA^f + e\ + m^] (t)\x) = = 1, 2, 3, 4), 

where Aj is the eigenvalue of F^^, corresponding to the eigenvector Fj. 

In the case of the coUinear electromagnetic field E = (0, 0, -E), B = (0, 0, 5), 



'-ri^'F., = iE 



^3 
^3 



B 



^3 
as 



and we can choose its eigenvectors and corresponding eigenvalues as follows: 



/ 1 \ 


1 



/ \ 

1 




F, 



Xi=iE-B X2 = iE + B A3 



/ 1 \ 



-1 

V / 

iE-B A4 



/0\ 

1 


V 1 / 

iE + B. 



(B.6) 



(B.7) 



(B.8) 
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B.l Free field 



First, we apply the method explained above to the free field case [A'^ = 0) and derive 

0. Because 



solutions of the Dirac equation 



*14 



In this case, we can choose any independent vectors as Fj since i^l'^l'^F^y 

2 = 1,2,3,4) 



the eigenvalue Aj is for any Fj, Eq.(B.6) becomes 

{dl + m^)(P\x) = Q 



(B.9) 



and its solutions are 

. . 1 



1 



[X) 



v/(27r)32^p V(27r)32^, 



giwpi+ip-x 



1,2,3,4). (B.IO) 



Define ^i,{x) as spinors made by Eq.(B.3) from the solutions (B.IO): 



\i0 + m 



(X)F,; 



[±70cUp-7-p + m] ±0^(x)F,. 



(B.ll) 



Notice that the operator in r.h.s is proportional to the projection operator to positive 
or negative energy state A-i-(p), respectively; [±7°ct;p — 7 ■ p + m] = 2mA-|-(±p). For 
any choice of Fj, the spinor of positive [negative] energy state is gotten from the pos- 
itive [negative] frequency solution of the "Klein-Gordon" equation (B.9). Choice of Fj 
affects only a spin state. 

To calculate the solutions explicitly, we use Fj defined by (B.8). It is sufficient to use 
only Fi and F2 because only four spinors are independent among eight spinors which are 
derived from four kinds of Fj. Spinors made from F3 and F4 can be expressed by linear 
combination of those from Fi and F2. 

Substituting the exphcit form of Fi and F2 into Eq.(B.ll), and rename ^^{x) as 
^pf{x) and ±'?/'p(a;) as ±'?/'p|'^(x), then we obtain 
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; , F2 + Juj.p - p^Fa 

V^p - Pz V^p " Pz 



1 



m 



m 

.V^p-P^ 
m 



■T2 + 



:F2 + 



Px + iPy 
y/ujp + Pz 

Px - iPy 

V^p - Pz 

Px - iPy 
^UJp + Pz 



^2 - a/^^p + Pz^Z 

ri + ^J^p - p^F. 



v/(2vr)32a;p 
1 

J v/(27r)32a7p 
1 



g-ia;pt+jp-x 



Fl - ^J^p^PzXi 



1 

V(27r)32cu/ 



g-iojpi+ip-x 



hiojpi+jp-x 



(B.12) 



*^^In the free field case, it is more general and convenient to derive solutions of the Dirac equation by 
boosting solutions in a rest frame. However, in the presence of an electromagnetic field, the boosting 
method is not convenient because an electromagnetic field is transformed by boost, and free field solutions 
derived by the method which is identical to that in electromagnetic fields makes the calculations in 



Sec 



clear. 
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Because these solutions satisfy the orthonormal conditions 

uv'iT,+<7) = (-<r,-C'^) = 5s.'5'(p^ (#^T,-<:') = o {s,s'=],i), (B.13) 

where the inner product is defined by 



(B.14) 



the quantized free field operator can be expanded by these solutions: 



(B.15) 



B.2 In electromagnetic fields 

We solve the Dirac equation in the presence of the collinear electromagnetic field E = 
(0,0,E),B = (0,0,5). Choice of gauge as A'^ = {0, -By,0, -Et) converts the "Klein- 



Gordon" equation (B.6) to 



[d^ - (di + leByf - dl - {83 + ieEtf + + eX^] (j)\x) 
This can be solved in the same manner with Sec J2.1.2i A solution is 



0. 



(B.16) 



KM 




Dniv) — e 



(B.17) 



where the subscript p is the abbreviated expression for {px,Pz,n) {n is the Landau level) 
and variables ^ and t] are introduced as 



-,{eEt + pz), r] 



eE' ' V eB 

The time-dependent part 0p(t) satisfies the equation 



(eBy + px). 



(B.18) 



d^^ 4 2eE 



<P'{t) = 0, 



(B.19) 



where Cn = (2?t- + l)eB. We choose two sets of solutions of Eq.(B.19) with the same 



motivation as explained in Sec, 2.1 



and 



{2eE)^ 
e ^ " 



(2eS)3 



5 —ia„~i- 



(t) 



{2eE)i 
e ^ " 



tD_ 



{2eE)i 



{2eE)i 



^out,l 



^out,2 



(t) 



{2eE)i 



(2eE)i 



tD* 



I. -L 



(B.20) 



(B.21) 
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where parameter a]^ and are defined by 



a. 



2eE 2eE _ n 1 . . . 



I _ rre^+Cn-e-B _ m^+{2n+2)eB 



(B.22) 



From these solutions, we can construct spinors using Eq.(B.3). Define [x) as 

which is derived by substituting Z{x) = N j^'^""^^'^ {t)Dn{ri)e'-''^'''^^'^^^^T i into Eq.(|B.3[) and 



±'^pi°"*(^) '^(^) ^ A^i(/)p''°^'''^()f:)-D„(r7)e**^^^^+P^^^r2 (A^ is a normalization constant). 
The indices t and | distinguish the spin state: t means spin parallel to the magnetic field 



and I means anti-parallel, which are known from Eq.(|B.22|). The resultant spinors are 



e 4' 




2VeE V 27r V vr y 



2^^ 




L /eBy 1 



2VeE V 27r V TT y 
/2eE 



'ail 



2^ 



g-f^n / I f eB\ 4 1 




2Ve^ \ 271 \ n J 



2^^ 




2VeE \ 2tt \ 7T J 



2^^ 

(B.23) 



and 




L feBy 1 



2 y/eE V 27r V TT y 
/2eS 



m 



D_Men)Dn{v)^i 



ar, 



e"4"" / L f eB\ * 1 
2v^ V 27r V TT y 



-V2eBD* , (e^^0^^n-i(^7)r2 - v^e^^Z^* (e^*0^n(^)r3 



3i(Pj;X+P22) 

2^r 



e 4 



feB\~^ 1 




^^-.ai(e^^0^n(r/)r2 



2 yeE V 27r V TT 



27r 



3i(Pa;a;+Pz2) 



2-K 
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L feBY 1 



2VeE V 27r V TT y 



These solutions satisfy the orthonormal condition 

(+V^ps, +^qs') = (-^ps, -^qs') = ^{Vx ' qx)S{p, - Qz) t^n,n,Sss' /aS = in, OUt 



oiiPxX+PzZ) 



2n 



(B.24) 



. (B.25) 



The quantized field operator in the presence of the coUinear electromagnetic field can be 
expanded in two ways: 



/ '^^^ J 




1 dp^ ^ 


[ dpz 



(B.26) 



We can obtain solutions in the pure electric field, = (0, 0, 0, —Et), from Eq.(B.23) 



and (B.24) by the replacements 



V2v 



(f ) ' ^V2^D^+M ^ {p, - zpy) 
(f ) ^ -}^V2^nD,,_,{r^) ^ (p, + tpy) 



2n 



an — > a — 2g£;- 



The results are 



2VeE 



+%/2^e^*Dl„(-e5*Or3 



gip X 



2VeE 



7(2^ 

jpx 

V(27r)3 



oipX 



e 4 



2Ve^ 

-mTeJ^D_,,_i(-e^^e)r4] 



2VeE L 



oipX 



(B.27) 



(B.28) 
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and 



2VeE 



e-4' 



2VeE 



e-4' 



0«PX 



2VeE 



-Cf (^) = 



o«pX 



y(2^- 



(B.29) 



Solutions in the pure magnetic field, = (0, — -By, 0, 0), can be obtained from the 



free field solutions (B.12) by the replacements 



/2tx 



1 



/' eS ^ i 1 

27r V 2-iT 



(B.30) 



v/m2+p2 + (2n + 2)e5 



The results are 



1 rZT/e^y 1 



=L'„(r7)ri - 



ujI -pz 



-Pz 



znDn-i{r])T2 



(27r)22cjj5 



1 [L /eBY 1 



v/2 V 27r V TT / 



ujI+Pz 



(27r)22cuj5 



ivl+Pz 



f^pi(^) 




1 / L {eB\ 4 1 



V 27r V TT y 



:I)„(r7)r2 + I^n+i(r?)ri 



(27r)22a;^ 



jLoi^t+i(p3;X+pz2;) 
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1 



,B / X 1 L feBY* 1 



uji+Pz y^^p+Pz (B 31' 



1 



(27r)22^J, 



Notice that, in this case, distinction between in and out solutions is unnecessary since 
the magnetic field affects only transverse part of the solution and time-dependent part 
g=FiLjpt unchanged, which is separable into positive and negative frequency in all time. 
In other words, a pure magnetic field cannot cause pair creation. 
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